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Summary of Work Completed 


The work completed to date is comprised of the following: a simple vehicle model 
representative of the aerospace plane concept in the hypersonic flight regime, fuel-optimal 
climb profiles for the unconstrained and dynamic pressure constrained cases generated 
using a reduced order dynamic model, an analytic switching condition for transition to 
rocket powered flight as orbital velocity is approached, simple feedback guidance laws for 
both the unconstrained and dynamic pressure constrained cases derived via singular 
perturbation theory and a nonlinear transformation technique, and numerical simulation 
results for ascent to orbit in the dynamic pressure constrained case. 

A hypersonic research airplane concept studied by NASA in the mid- 1970’s has been 
selected as a nominal vehicle configuration. A variety of windtunnel data is available for 
this vehicle in the open literature over the Mach range 0.2 to 6.0. A full-scale vehicle of 150 
feet total length and 200,000 pounds gross take-off weight representative of the X-30 
research aircraft is assumed. The vehicle model consists of separate modules for the 
estimation of aerodynamic and propulsive forces. The aerodynamic model, currently 
restricted to the hypersonic regime, is based on coefficients obtained by curve fits to 
windtunnel data at Mach 6.0. A parabolic drag polar is assumed and the induced drag 
coefficient is assumed independent of Mach number via the Mach Independence Principle. 
A largely analytic aerodynamic model based on a combination of Newtonian flow and blast 
wave theories was also investigated. The propulsion system is assumed to consist of a 
bank of six SCRAMJET engine modules that operate continuously at stoichiometric 
conditions when above Mach 10 and a rocket engine rated at 15,000 pounds gross thrust 
that can be turned on or off as dictated by optimality conditions. A simple conceptual 
SCRAMJET model, largely analytic and well suited for trajectory optimization, has been 
developed and is outlined in this report. 

The mission of single-stage ascent to orbit was considered, and fuel-optimal ascent 
trajectories were generated numerically using a reduced order dynamic model. This model 
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results from the following: application of an energy-state approximation, the assumption of 
a non-rotating Earth, the assumption that flight is constrained to a vertical plane, and finally 
by imposing time scale separation of altitude and flight path angle dynamics from energy 
and mass dynamics via singular perturbation theory. Rocket throttle control as currently 
modeled appears linearly in the Hamiltonian. The possibility of a singular arc in rocket 
throttle setting was examined and was shown to be non-optimal. Thus, the optimal control 
for rocket propulsion is bang-bang. An analytic condition for switching rocket thrust on, as 
orbital velocity is approached, was derived and is presented in this report A zeroth-order 
boundary layer solution which accounts for altitude and flight path angle dynamics was 
formulated. A simple costate approximation was adopted to allow implementation but 
resulted in an unstable feedback guidance law. An improved costate estimate was 
subsequently formed by linearizing the boundary layer necessary conditions about the 
reduced solution. The unstable modes were then successfully suppressed. The resulting 
stable nonlinear feedback guidance law applies to the unconstrained case. The incorporation 
of a dynamic pressure constraint leads to difficulties in generating a boundary layer 
solution. A suboptimal lift control solution was derived via a nonlinear transformation 
technique as an alternative approach while the afore mentioned case is being investigated. 

A numerical simulation of the hypersonic phase of fuel-optimal ascent-to-orbit was 
carried out using the suboptimal guidance law mentioned above for the dynamic pressure 
constrained case and is documented in Section 6 of this report. Lastly, future research 
objectives are recommended. These include: theoretical investigation of the state 
constrained boundary layer problem, extension of the vehicle model to the subsonic and 
supersonic flight regimes, the modeling of SCRAMJET thrust dependence on angle of 
attack, modeling of the component of thrust that contributes to vehicle lift, and 
investigation of optimality conditions associated with discrete variations in SCRAMJET 
thrust. Also recommended is the consideration of additional constraints (temperature, 
acceleration and lift limits) and examination of three-dimensional maneuvers as may be 
required for plane change, lift modulation, or mission abort. 
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1. Introduction 


The Space Transportation Problem 

Recent studies [ 1 ,2] which review the state of space transportation warn that a cheaper, 
more reliable means for transporting both people and cargo to and from earth orbit must be 
developed in the next 20 years if the United States is to maintain its position as a world 
leader amongst space-faring nations. The current U.S. Space Transportation System 
(Space Shuttle Fleet) represents an effort to build one vehicle to serve many roles. Despite 
its technical success, it is unlikely that a future derivative of this vehicle can achieve the 
operational efficiency required to remain competitive in the growing international space 
launch business or enable the U.S. to open the so called "space frontier". Studies indicate 
that our projected transportation needs will best be served by a mix of expendable and 
reusable vehicles. Specifically the functions of cargo transport to orbit and two-way 
passenger transport should be separated. In the case of either mission we require fully 
reusable vehicles that are robust and reliable. Numerous configurations, fuels, propulsion 
methods, launch modes, and other characteristics have been studied to determine the more 
promising approaches. Two types of vehicles have emerged: vertically launched rocket 
vehicles of both one and two-stages and horizontal-take-off single-stage-to-orbit air- 
breathing vehicles. The latter are commonly referred to as aerospace planes [3]. These two 
competitors, advanced rocket and aerospace plane technologies, both promise drastic 
reduction in the cost of achieving orbit and must be developed now if we are to be ready to 
exploit the opportunities of space in the 21 st century. In this report, the problems 
associated with the guidance and control of the aerospace plane configuration are examined 
in light of its desired operational objectives. Then, the progress made in deriving guidance 
laws suitable for implementation on-board such a vehicle is reported. 

The Promise of NASP Technology 

The aerospace plane concept involves winged or all-body vehicles, fueled by liquid 
hydrogen, that can depart and land horizontally from conventional jet runways. The 
configuration of principal interest would be capable of flying to low-earth orbit using only a 
single stage. The critical technologies that must be advanced include air-breathing 
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supersonic combustion ramjet (SCRAMJET) engines, high temperature materials, and 
hypersonic configuration aerodynamics [1,4]. This technology, if developed, would aid 
not only transportation to low earth orbit, but also a host of other potential hypersonic 
missions, both military and civilian, identified in various government sponsored studies 
over the past several decades. Configurations could range from advanced interceptors and 
high performance reconnaissance aircraft to transports capable of cruise at Mach 6-12[4-7]. 

A successful aerospace plane, configured for ascent to low-Earth-orbit, would not only 
achieve a payload mass fraction an order of magnitude greater than current rocket systems, 
but would also do the following: eliminate conventional rocket staging, offer on-demand 
assured launch for space rescue and national security, provide greater basing flexibility 
through operations from conventional airfields and self-ferry flight, and improve system 
survivability by eliminating our reliance on just two U.S. launch complexes. Such a vehicle 
could also provide improved mission safety through multi-engine redundancy, aircraft-type 
control, abort capability, and alternate mission-recovery paths (cross-range flight 
capability). Airplane-like operations would greatly reduce the large number of ground 
support personnel and eliminate vertical-assembly buildings, launch pads, special recovery 
flight operations, solid boosters, external tanks, and other logistics burdens that 
characterize the Shuttle. In short aerospace plane technology, when mature, could offer 
efficient, reliable and economic access to orbit [8]. 

The History of the NASP Program and Its Program Objectives 

The aerospace plane concept, which dates back to the 1950’s, was seriously 
investigated in the United States during the 1960’s. Development was abandoned at that 
time due to technical barriers. Subsequently, pure rocket technologies began to dominate 
our research efforts. All but a few research programs in SCRAMJET propulsion at NASA 
Langley and with the Navy were terminated. In 1982 the Defense Advanced Research 
Projects Agency (DARPA) again began investigating the limits of air-breathing propulsion, 
aerodynamics, materials, and structures in an initiative termed Copper Canyon. In the 
period between 1982 and 1985, DARPA redefined the aerospace plane concept with 
laboratory support from NASA, the Air Force, and the Navy. By 1985 it was determined 
that advances in aerodynamics, structures and air-breathing hypersonic propulsion had 
significantly lowered the technological barriers encountered earlier, and a decision was 
made to initiate a technology development program [4,7,9]. This program is referred to as 
the National Aerospace Plane (NASP) Program, and features wide participation by 
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Government agencies and industry. It targets the maturation of key technologies, and plans 
the feasibility demonstration of a radically advanced engine. The preliminary design of the 
airframe needed for an experimental flight research vehicle analogous to past X-l and X 15 
research aircraft is underway. The new engine will be built and ground-tested up to about 
Mach 8, while the new structures and materials needed to fabricate such a vehicle are being 
developed and tested. Since no ground test capability is available to simulate full scale flight 
conditions at Mach numbers much above this, and indeed none seems feasible, the 
propulsion system and its integration into the airframe must be validated experimentally in 
flight. A research aircraft, dubbed the X-30, is scheduled to fly in 1994-95 in a test 
program that will demonstrate hypersonic cruise and acceleration into low-Earth orbit [1,4]. 
If the flight program confirms the feasibility of aerospace plane technology, definition and 
full-scale engineering development of next generation transport vehicles can be initiated 

using this vehicle concept 


Guidance, Navigation and Control of the Aerospace Plane 

This research project addresses the problems associated with the guidance and control of 
vehicles that may evolve from the NASP program. These vehicles, although varied in size 
and detail according to the intended mission, will be of similar character where guidance, 
navigation and control (GN&C) issues are concerned. This family of vehicles will be 
referred to as "the" aerospace plane throughout this report 

The GN&C of the aerospace plane will ultimately include many aspects, including 
ascent to orbit aero-assisted maneuvers, and the like. An essential component to achieving 
the transportation cost reductions mentioned earlier is that of drastically simplifying launch 
operations. In order to make possible an order of magnitude reduction in the cost per 
pound of payload placed in orbit, aerospace plane operations will have to approach those of 
current commercial airlines [1] . Current GN&C technology, with its reliance on pre- 
mission planning is inadequate to meet the challenge of automatic and adaptive trajectory 
control [9,10]. On-board guidance algorithms are needed to provide rapid airline-like 
operations and to respond to changing mission demands. This research effort seeks to 
develop computationally efficient and robust analytical and computer methods suitable for 
on-board flight trajectory optimization. Energy methods and singular perturbation theory. 
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which have been successfully applied to similar problems for fighter aircraft [1 1-14], are 
the principal tools to be applied in this endeavor. 

Review of the Document that Follows 

What follows is a report of the progress made toward this end during the period July 1, 
1987 to October 31, 1988. Section 2 of this report contains a description of the selected 
vehicle configuration and its intended mission. A hypersonic aerodynamic model and 
models for SCRAMJET and rocket propulsion systems are presented in Section 3. Section 
4 addresses the issues of vehicle sizing and the selection of engine scaling parameters. 
Section 5 provides an introduction to the ascent-to-orbit guidance problem and, after a 
suitable set of dynamic equations are introduced, various assumptions which lead to a 
reduced order model are investigated. Next, singular perturbation theory is used to derive a 
means for approximating the minimum fuel climb path and to develop guidance laws in 
feedback form. A nonlinear transformation technique is also employed as an alternative to 
solving the boundary layer problem. The possibility of intermediate values of rocket throttle 
exists and an argument is presented to show that such settings are non-optimal. An analytic 
form for the rocket throttle switching condition is also presented. The approximate fuel- 
optimal climb paths for the unconstrained and dynamic pressure constrained cases are 
presented in the text as the analysis is developed. A numerical simulation of the hypersonic 
phase of ascent-to-orbit using a derived guidance laws for the dynamic pressure 
constrained case is documented in Section 6. The conclusions for this effort and 
recommendations for future research are given in Section 7. Appendices, which further 
detail various technical issues, appear at the end of the report. 
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2. Vehicle and Mission Description 


The aerospace plane concept in general involves winged vehicles, fueled by liquid 
hydrogen, that can depart and land horizontally from conventional runways, maintain 
hypersonic cruise in the upper atmosphere for long durations and accelerate to orbital 
velocity. All-body aircraft have also been considered as candidate vehicles. It is not clear at 
this time which vehicle type will be most suitable. The high ignition speed of hydrogen fuel 
is required for supersonic combustion and this fuel’s greater energy content is 
advantageous, but it also has a higher specific volume than conventional hydrocarbon 
fuels. The resulting high fuel volume requirement may dictate a blended wing-body or an 
all-body design. Since these issues are yet to be resolved, this research effort will assume a 
winged vehicle. Many potential missions have been identified for such vehicles [6]. The 
particular mission of single-stage ascent-to-orbit which promises, by the use of air- 
breathing hypersonic propulsion and greatly reduced launch operations, order of magnitude 
reductions in the cost of placing payloads in low-Earth orbit is especially attractive. The 
National Aerospace Plane Program seeks to develop a hypersonic research aircraft, 
designated the X-30, which is capable of demonstrating the feasibility of the technology 
required to achieve such operations. The research aircraft will be designed to realize four 
primary goals during flight tests [15]: 

• To simulate cruise for extended durations at speeds between Mach 5 and 
Mach 10 and altitudes well above 100,000 feet 

• To demonstrate operation into and out of ordinary airports, including 
environmental acceptability. 

• To demonstrate flight into orbit from a runway, powered by air-breathing 
engines and carrying several thousand pounds of payload. 

• To establish rapid turnaround in space operations. (Considered by some 
even more important than maximizing the payload mass fraction to orbit.) 
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Clearly the X-30 , if successful, shall possess most if not all of the aerodynamic and 
propulsive characteristics of the general aerospace plane concept. Thus we shall consider a 
vehicle model representative of the X-30 as suitable for investigating problems associated 
with the GN&C of this family of vehicles. Unfortunately, the X-30 development program, 
currently in the technology demonstration phase, has been classified as secret by the U.S. 
Department of Defense. As a result, few details will be forthcoming in the open literature. 
This is not, however, a serious obstacle, since a number of hypersonic research vehicles of 
similar nature were designed in the 1970’s and, accompanied by suitable windtunnel data, 
appear in the open literature. 

A vehicle representative of the X-30 has been selected for this project from the series of 
hypersonic vehicle configurations studied by NASA in that time period. This aircraft, for 
which three-dimensional windtunnel data are available over the Mach range 0.2 to 6.0 [16- 
19], is a hypersonic research vehicle concept with a 70 degree swept delta planform. It 
features an airframe integrated SCRAMJET propulsion system. It was designed to be air- 
launched from beneath the wing of a B-52 much the same as the earlier X-15. A 
photograph of a 0.021 scale model of this vehicle and the configuration’s general 
dimensions, normalized to body length, are reproduced from reference [18] as Figures 1 
and 2, respectively. Various configurations of this particular vehicle were tested. The 
configuration designated B!W lf V t F d E was selected as the most appropriate for this study. 

Table 1 presents a summary of the selected configuration features. 


Table 1 Summary of Vehicle Configuration Features 


B i Body one, with a high profile nose reflecting a forward cockpit location 

Wjf Wing one, 70° swept delta planform with positive camber, forward location 
V t Vertical wing tip fins (as opposed to a center vertical tail), with 7.5° toe-in 
F d Additional forward delta wing 

E Underslung SCRAMJET engine modules 



ORIGINAL PAGE IS 
OF POOR QUALITY 



Figure 1: Configuration Model. 





11 


This design, based on a fixed geometry modularized SCRAMJET concept that is 
integrated closely with the basic airframe, is typical of the hypersonic air-breathing systems 
studied over the past two decades. The multiple engine modules are attached to a forebody 
precompression surface and exhaust over an aft body-nozzle surface. The inlets of the 
multiple rectangular SCRAMJET modules efficiently capture precompressed airflow 
contained between the vehicle and the forebody shock wave. The afterbody nozzle serves 
to increase SCRAMJET nozzle expansion area and allows the external SCRAMJET nacelle 
to be nearly stream aligned at the design Mach number for maximum installed thrust 

performance [20]. 

SCRAMJET engines, however, are incapable of operation below Mach 3 or 4 and 
suffer greatly reduced thrust capability at extreme altitudes. Aerothermodynamic constraints 
and the desire to achieve orbit will eventually force the vehicle to climb to such altitudes. 
Rocket propulsion is required for attitude control in space and re-entry. For these reasons 
the vehicle design must incorporate a multi-mode propulsion system. Candidate engine 
types include turbojets, ramjets, air-turbo-ramjets, SCRAMJETs, and rocket engines. One 
NASP Program objective is to achieve orbit using air-breathing propulsion alone, but many 
doubt the SCRAMJETs ability to power the vehicle all the way to Mach 25 [21]. Thus it 
is of interest to consider the use of a rocket engine prior to exiting the Earth's atmosphere. 
Which combination of engines is optimal and when to transition from one type to another 
are key questions to be answered. This report is focused on the hypersonic flight regime 
and addresses the question of SCRAMJET / rocket transition in a later section. 

The above configuration, scaled to a length of 150 feet and weighing 200,000 pounds 
when fully fueled with liquid hydrogen, is used throughout the remainder of this report as a 
vehicle model representative of the X-30 research aircraft [22]. The resulting reference area 
used for defining aerodynamic coefficients, (taken to be the projected area of the wing 
planform, including the part encompassed by the body) is 3780.0 square feet. In the next 
section, consideration is restricted to the hypersonic flight regime and a dual-mode 
propulsion system consisting of SCRAMJET and rocket engines is defined. Modeling of 
the vehicle's hypersonic aerodynamic characteristics and of the dual-mode propulsion 
system is discussed in detail. More about the sizing of the vehicle and its propulsion system 
is given in the section that follows. 
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3. Vehicle Modeling 

As detailed in the previous section, a hypersonic research airplane concept studied by 
NASA in the mid 1970's was selected as a nominal vehicle configuration. This design is 
based on a fixed geometry modularized SCRAMJET propulsion system that is closely 
integrated with the airframe. Consideration has been restricted to flight above Mach 5 and a 
dual-mode propulsion system (SCRAMJET / rocket) has been assumed. Vehicle sizing, 
discussed in the next section, results in a full scale vehicle of 150 feet total length and 
200,000 lbs. gross take-off weight. For the purpose of trajectory optimization and 
simulation of derived guidance laws, a model reflecting the aerodynamic and propulsive 
characteristics of the vehicle is required. This model must be of a complexity commensurate 
with the expected accuracy of the performance analysis. In determining the nature of 
optimal flight profiles it is necessary that the models correctly predict the qualitative 
behavior of the aircraft, although to a limited degree, quantitative accuracy can be 
sacrificed for the sake of simplicity. In the study of hypersonic aircraft, a need currently 
exists for simple performance codes that use (to the extent possible) analytic aerodynamic 
and propulsion models and efficient optimization algorithms [7]. In this section the 
aerodynamic and propulsive models used to generate numerical results are presented. 
Included is the rationale for the models chosen and an indication of their shortcomings. 
Appendix A reviews the progress made in developing an analytic aerodynamic model for 
slender-bodied aircraft in the hypersonic flight regime. 

3.1 Aerodynamic Model 

A variety of windtunnel data are available for the selected vehicle configuration in the 
open literature over the Mach range 0.2 to 6.0 [19]. The aerodynamic reference area (s), 
taken to be the projected area of the wing planform, including the part encompassed by the 
body, is 3780 sq. ft. The angle of attack along the optimal climb path will be shown to 
remain quite small, thus it is reasonable to assume a parabolic lift-drag polar of the form 


C d = C do + KCj 2 


( 1 ) 
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atmosphere and/or operating as space launch vehicles for delivering payloads to orbit. 
Supersonic combustion ramjet (SCRAMJET) engine technology is the key to making such 
a vehicle concept a reality and when closely integrated with a host of various other 
advanced technologies may make possible a fully-reusable single-stage-to-orbit launch 
vehicle with a payload mass fraction an order of magnitude greater than that of current 
rocket systems [8]. 


Overall engine performance can be characterized by average specific impulse (Isp) 
which is defined as the number of units of thrust produced per unit of fuel weight flow 
rate. The units of Isp are seconds and the larger the value the more efficient the propulsion 
system. Rockets are limited typically to less than 500 seconds Isp, while a multi-mode 
aerospace plane propulsion system incorporating SCRAMJET engines is expected to 
average between 1,500 and 2,000 seconds Isp [8]. The potential performance gains are 
indicated in terms of specific impulse in Figure 5, which was reproduced from reference 
[6]. Clearly the ramjet and supersonic combustion ramjet (SCRAMJET), both airbreathers, 
can provide (as seen in Figure 5) efficient cruise propulsion for hypersonic vehicles. Note 
that the use of hydrogen rather than hydrocarbon (e.g. kerosene) fuels can improve engine 
performance at all flight speeds (see Figure 5). In particular, the liquid-hydrogen-fueled 
SCRAMJET offers the potential of Mach 7 performance comparable with that of 
hydrocarbon supersonic turbojet [6]. 


3.2.1 SCRAMJET Model 

The supersonic combustion ramjet (SCRAMJET) has been well understood 
conceptually since the early 1960's but has been experimentally validated only in ground 
tests, and only up to Mach numbers of about 8. To be attractive for orbital launch, its Mach 
number range must extend at least to 12, and preferably beyond 15 [1]. The lack of 
appropriate unclassified experimental data, cycle analyses, and combustion analyses 
requires the use of a simple conceptual model for the purpose of vehicle trajectory 
optimization. What follows is a brief description of the model being used in this research 
effort and the philosophy behind it. 



RAMJETS 

(VARIABLE 
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Figure 5. Propulsion Options. 
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Conceptually, the SCRAMJET is as simple an airbreathing combustion device as one 
could imagine. In the case of the vehicle concept oudined in Section 2, the entire underside 
of the vehicle plays a role in the operation of the propulsion system. Figure 6 shows the 
basic configuration. Mechanically, the device can be thought of in terms of three elements. 

These are: 


1. Diffuser 

2. Combustor 

3. Expansion nozzle 

Hypersonic vehicle designers attempt to utilize the forward fuselage, strakes, and wings to 
provide the majority of the diffusion. The lower part of the three-dimensional oblique 
shock formed at the leading edges is tailored to the shape of the combustor inlet so that air 
enters at approximately Mach 3, but this depends on the flight speed. Combustion of 
hydrogen fuel takes place in the duct at supersonic speeds in order to minimize energy 
losses due to dissociation, which would be enormous if the more conventional subsonic 
ramjet cycle were to be used in high speed flight. Liquid hydrogen is the fuel of choice not 
only because of its high energy content, but because it can be made to bum in a supersonic 
flow due to its wide flammability limits and high flame speed. Finally, the combustion 
products are expanded through a nozzle, which, like the diffuser is designed into the 
contour of the lower fuselage. 

As Figure 6 suggests, the propulsion system is mosdy diffuser and nozzle. While these 
elements are fairly easy to model from the thermodynamic cycle point of view, the 
aerodynamics are quite complex, giving rise to a challenging design problem. 
Computational fluid dynamics (CFD) numerical techniques are being relied upon in 
conjunction with a new family of hypersonic test facilities to yield practical design 
solutions. Unfortunately, information on the current research is classified, so that realistic 
design data is not available for projects of this type. 

The computational model used here to represent the SCRAMJET propulsion system was 
deliberately designed to be readily updated as new information becomes available. It 
directly accesses a standard atmosphere model (also easily adjustable to provide non- 
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Figure 6. Schematic Layout of SCRAMJE r Propu 
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standard opera, ing condi, tons), which simplifies its incorporation into a trajectory 
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standard thermodynamic models or by means of optimal destgn curve fits such 
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performance and should be properly modeled. 

Figures 7-10 present calculated SCRAMJET perfotmance as a function of Mach number 
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“r ratio to remain a, or below its stoichiometric value. The resulting change ,n slope 
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in the thrust curve that occurs as f reaches its stoichiometric value will later be shown to 
influence the unconstrained fuel-optimal climb path. The SCRAMJET computational 
algorithm also provides estimates of fuel flow rate, propulsion module thrust, propulsion 
drag, and the standard performance parameters. Additional references on SCRAMJET 
propulsion systems are included in the bibliography as [28-35]. 

3.2.2 Rocket Model 

In addition to its SCRAMJET engines, the X-30 must carry rocket propulsion for 
attitude control in space and reentry. It is also of interest to consider the use of a rocket 
engine prior to exiting the Earth's atmosphere. The Pratt & Witney RL10 rocket engine 
used on the Centuar upper stage is deemed suitable for this purpose [21,36]. The RL10 is 
rated at 15,000 lbf thrust at 200,000 ft. with a nominal specific impulse of 444 sec. 
Propellants are liquid oxygen and liquid hydrogen with a nominal oxidizer-to-fuel ratio of 
5:1. This engine, capable of multiple starts, has been tested in advanced versions that 
include variable thrust. The rocket model used in this study assumes the performance of the 
RL10 and makes a simple thrust correction below 200,000 ft for losses due to atmospheric 
back pressure [37], Thrust of the rocket is given by 

T r = T vacum — A e p(h) (2) 

where A E represents nozzle exit area and p is atmospheric pressure given as a function of 
altitude. Figures 1 1 through 14 present predicted rocket performace as a function of altitude 
in terms of fuel specific impulse, thrust, thrust specific fuel consumption, and mass flow 
rate given a nozzle exit area of 1 .0 square foot. 
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Altitude in feet 

Figure 1 1. Predicted rocket performance, fuel specific impulse 
versus altitude, for a nozzle exit area of 1.0 square foot. 



Altitude in feet 

Figure 12. Predicted rocket performance, thnist versus altitude, 
for a nozzle exit area of 1.0 square foot. 


TSFC (slug/ hr/ 1 bf) 
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Altitude in feet 


Fieure 13. Predicted rocket performance, thrust specific 
fuel consumption versus altitude, for a nozzle exit area of 
1.0 square feet 



Figure 14 Predicted rocket performance, mass flow rate 
versus altitude, for a nozzle exit area of 1.0 square feet 
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4. Vehicle Sizing 


As previously detailed, a hypersonic research airplane concept studied by NASA in 
the mid-1970’s was selected as a nominal vehicle configuration. This design is based on a 
fixed geometry modularized SCRAMJET propulsion system that is closely integrated with 
the airframe. Consideration has been restricted to flight above Mach 5 and a dual mode 
propulsion system (SCRAMJET /rocket) has been assumed. This vehicle has been sized to 
resemble the scale of the proposed X-30 research aircraft and results in a full scale vehicle 
of 150 feet total length and 200,000 lbs. gross take-off weight [38]. The resulting 
aerodynamic reference area (s), taken to be the projected area of the wing planform, 
including the part encompassed by the body, is 3780 square feet. 

The hypersonic aerodynamic characteristics of this vehicle configuration were used to 
size the SCRAMJET engine inlet area while assuming no other type of propulsion system 
to be operating. Figure 15 presents the vehicle's calculated level flight envelope as a 
function of SCRAMJET engine inlet area assuming weight remains constant at the assumed 
take-off value of 200,000 pounds and that lift equals weight. The inner-most envelope 
corresponds to 360 square feet of inlet area (6 SCRAMJET modules at 60 square feet each) 
which provides sufficient thrust for cruise between Mach 8 and Mach 12. The inlet area 
required to cruise increases dramatically with Mach Number, approaching the requirement 
for 900 square feet as orbital velocity is approached. This case corresponds to the outer- 
most envelope shown. This trend agrees with that indicated in reference [39]. Note that a 
maximum allowable dynamic pressure constraint (q=2000 psf is shown in the figure) 
severely limits the accessible flight envelope and thus greatly constrains the ability to 
optimize the ascent trajectory. The inlet area required per engine module exceeds 100 sq. 
ft. if the level flight envelope is to cover the Mach range from 5 to 20. This value is 
approximately 10 times the inlet area (note that inlet area refers to the cross sectional area of 
the combutor inlet) of the configuration shown in Fig. 2. Fuselage drag predictions based 
on the referenced windtunnel data have not been modified to reflect this change in vehicle 
frontal area and it has been assumed that the ram drag calculations made in computing 
SCRAMJET thrust account for the additional drag incurred. A more accurate analysis 



26 


would necessarily have to correct the fuselage drag estimate for this change in 
configuration. 


Figure 16 again presents the calculated level flight envelope as a function of inlet area 
given only SCRAMJET propulsion, but now it is assumed that lift equals weight minus 
centrifugal force. Weight is again held constant at its assumed take-off value of 200,000 
pounds. This figure indicates that for the postulated vehicle configuation a narrow corridor 
is avaiable for ascent to low-Earth orbit using SCRAMJET propulsion alone. Figure 17 
presents the equivalent SCRAMJET envelope but for weight assumed constant at 77,000 
pounds. This weight corresponds to an estimate of the remaining vehicle mass upon 
achieving orbital velocity. Having accounted for the decreased vehicle weight we find that 
the afore mentioned corridor is widened. 

Since rocket propulsion will be required for orbit circularization, on orbit maneuvers 
and initiation of reentry, it is of interest to examine the optimality of using this available 
rocket propulsion during ascent. Figure 18 presents the calculated level flight envelope for 
seperate SCRAMJET or rocket propulsion , and a combination of the two. The innermost 
envelope corresponds to a propulsion system consisting of 6 SCRAMJET engine modules 
with 100 square feet of projected inlet area per module. When compared to the curve 
corresponding to 6 SCRAMJET modules with 150 square feet of projected inlet area each, 
one sees that there is a level of total SCRAMJET inlet area (that is inlet area per module) 
between 100 and 150 below which orbit cannot be achieved. Increased engine inlet area 
comes with a tremendous drag penalty at the lower Mach numbers, thus it would appear 
that the use of rocket propulsion to agument SCRAMJET thrust would be advantageous to 
a point. A careful design trade-off in sizing the SCRAMJET and rocket engines, rather than 
simply sizing the rocket to perform its on-orbit duties, will produce the most efficient 
configuration. No such trade-off is attempted here. All of the numerical results that follow 
in this report were generated assuming 6 SCRAMJET modules with 150 square feet of 
projected inlet area per module and a rocket rated at 15,000 pounds thrust at 200,000 feet. 
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Figure 16. Level Right Envelope (SCRAMJET only) as a function of engine inlet area. 
Values shown are the inlet area of one of the six engines generating thrust. Again weight 
is constant at 200,000 lbs but now lift equals weight minus centrifugal force. 
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5. Trajectory Optimization and Guidance Law 

Development 


For the remainder of this report we shall consider only one of the many tasks that will 
be required of the aerospace plane's guidance and control system, that of ascent to low- 
earth orbit. Such missions are routinely performed by rocket powered expendable, or in the 
case of the space shuttle, partly reusable vehicles. The NASP technology program 
emphasizes the goal of full reusability because of its pronounced impact on cost and 
operational flexibility. In fact, analyses of aerospace plane operations and support costs 
have shown that the greatest cost savings come from reduced turnaround time [8]. 
Advanced guidance and control systems are a key component to establishing rapid airline- 
like launch operations. In order to reduce the cost of designing flight profiles, reduce the 
time required to respond to a changed payload or mission requirement and to improve the 
vehicle performance, on-board real-time optimal trajectory calculations will be required. 
The software must anticipate all possible mission requirements in order to avoid the need 
for mission dependent software modifications that lead to extensive preflight software 
testing requirements. Application of the necessary conditions for optimality in solving the 
trajectory optimization problem in general leads to a two-point boundary-value problem 
(TPBVP) that is difficult to solve. While some success in designing a reliable iterative 
algorithm to solve a TPBVP in an on-board computer has been achieved for orbit transfer 
[10], the diverse mission requirements of a general purpose aerospace plane will likely 
require that structured methods for order reduction be employed. 

Energy state approximations and singular perturbation methods have proven to be useful 
in deriving on-board trajectory optimization algorithms. These methods also contribute 
considerable insight into the nature of the optimal profiles and their relation to vehicle 
aerodynamic and propulsion characteristics. Most of the studies performed thus far have 
been devoted to fighter aircraft performance optimization in the context of minimum time 
intercept [ 14,40,41], These studies culminated in a series of piloted simulation evaluations 
at NASA Langley [42] and flight test demonstrations at NASA Dryden [43]. The 
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techniques that were used are currently being applied to optimal orbit transfer maneuvers in 
the upper atmosphere [44], 

Many of the modeling approximations employed for analysis of subsonic and 
supersonic aircraft optimal trajectories are not valid for a vehicle with hypersonic cruise and 
orbital capabilities. These can include the assumption of a flat Earth, constant vehicle 
weight, and a constant gravitational field. The problem of optimal airbreathing ascent to 
orbit is further complicated by the presence of severe trajectory constraints and the 
requirement for multiple modes of propulsion. SCRAMJET engines are incapable of 
operation below Mach 3 or 4. Alternate forms of propulsion are needed for take-off and 
landing, acceleration to supersonic speeds, and flight at the limits of the sensible 
atmosphere. The scope of this report is limited to flight above Mach 5; thus we shall not 
consider the required propulsion system transition in the supersonic regime. We shall, 
however, consider a transition in the hypersonic flight regime. As the vehicle accelerates 
toward orbital velocity it must continually gain altitude in order to avoid excessive 
aerodynamic heating. SCRAMJET performance degrades as altitude increases. The vehicle 
must convert to rocket propulsion as it leaves the atmosphere but the transition from 
SCRAMJET to rocket propulsion may well be advantageous prior to achieving orbital 
velocity. This transition to rocket propulsion during acceleration to orbital velocity must 
occur even if heating and other constraints are ignored since the fuel efficiency of a 
SCRAMJET degrades with increasing Mach number and will eventually fall below that of 
the rocket. The optimal point of propulsion system transition is identified in this report. 

5.1 Problem Formulation 

A good question to ask at this point is in what sense is the trajectory to be optimized? 
For the mission under consideration, namely ascent to orbit, one suitable goal is to 
minimize the total energy expended to achieve orbit. Since optimal space trajectories are 
well understood [45], we need only consider that portion of the flight within the sensible 
atmosphere, taken to be h < 259,000 feet. The minimum-fuel climb path is the selected 
goal for this analysis. 
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5.1.1 Dynamic Model 

Consider the equations of motion governing three-dimensional atmospheric flight of a 
point mass over a spherical rotating earth that are given below [46,47]. This set of 
equations is of an adequate complexity for our most detailed performance analysis. A 
stationary atmosphere and an inverse squared gravity field are assumed. 
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m = -f (r,V,rj) 


(9) 


The dots on the left-hand sides of equations (3-9) denote differentiation with respect to 
time. The state variables are: radius from the center of the Earth, r, longitude and latitude, 0 
and <(> respectively, the flight velocity, V, flight path angle and heading angle, Y and V 
respectively, and vehicle mass, m. The control variables are engine throttle, T\, vehicle lift, 
L, and bank angle, o. Figures 19 and 20 define the coordinate system and illustrate the 

geometric relationships between these variables. 

Lift and drag are defined as: 


L = Lift = l/2pV 2 Qs where Q = Q (a,Moo) (10) 

D = Drag = l/^pX^Qs where C d = C d (a,Moo) (11) 

The lift and drag coefficients, C, and C d , are assumed to exhibit known functional relations 
to the vehicle angle of attack, a, and the free stream Mach number, Moo- Furthermore, it is 

assumed that drag can be expressed in a conventional parabolic form as follows. 


D = qsC do + KL 2 /qs 


q = pV 2 /2 (12) 


where the symbol q represents dynamic pressure, s an aerodynamic reference area, C do the 
zero-lift drag coefficient, and K the coefficient of the induced drag component. Altitude (h) 
above mean sea level (r G ) is given by. 


h = r - r Q 


(13) 
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proceeds Other constraints which wili need to be incorporated into the analysis a. a later 
date include acceleration limits, lift limits, and aerodynamic heating tarns. A statable mode 
for aerodynamic heating has been identified in reference [49]. 


5.1.3 Simplifying Assumptions 

With suitable models defined for the prediction of aerodynamic and propulsive forces 

and the vehicle's dynamic response , i, is possible to simulate, using ^ ^ 

the vehicle's flight given a control history, or to seek a control h. story whtch wtll oputmze 
Tperflance index of interest. Rather than seeking solutions to trajectoty op- 
problems using the complex dynamic mcrfel described above we wish ^ 

of reduced order that will yield an acceptable approximate solemn. Thu ;* fte ■ ^ 
tha, the resulting nonlinear optimal control problem proves to be very d.ff.cult to dvo 
Both direct and indirect methods of solution are available, but each proves 
computationally intensive. Furihermore, none of these solution methods leads to a guidance 
law in feedback form. For these masons such methods are no, suitable for on-board rea - 
dme implementation and we must turn ,0 an approximate solution method. The energy s 
approximation has proven to be most useful in this regard (.1, ^ 

other simplifying assumptions normally employed for model order reduction when 
1 'tiering transport or fighter aircraft may no, be valid for a vehicle with hyperson.c 
capabilities. These assumptions typically include [46]: 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 
(7) 


Non-rotating Earth. 

Flat Earth. 

Constant acceleration due to gravity. 
Constant mass. 

Flight constrained to a vertical plane. 
Flight path angle small. 

Angle of attack small. 
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(8) Thrust aligned with the velocity vector. 

(9) Introduction of specific energy as a state variable 


By adopting equations (3-9) it has already been assumed, without justification, that we 
have: 

(10) Spherical Earth. 

(11) Stationary atmosphere. 

(12) Gravity field proportional to the inverse of the radius squared. 

(13) The vehicle may be modeled as a point mass. 

(14) Thrust vector constrained to lie in the plane of symmetry. 

(15) The location of the vehicle center of mass is fixed (affects trim and 
thus the calculation of drag and elevon settings). 


Let us examine each of the assumptions (1-9) individually , determine which may be 
employed for a hypersonic vehicle, and, where possible, determine the magnitude of the 
error introduced by doing so. Consider first those assumptions which lend themselves to 
analytical investigation, such as (1-3). 

Non-rotating Earth. Consider those terms in the given equations of motion that 
involve the rotational velocity of the earth, CD. This rotation gives rise to two forces, or 
accelerations. The first, known as Coriolis acceleration, gives rise to terms involving 2coV, 
which have an important impact on high speed, long range flight. The second, termed 
transport acceleration, gives rise to terms in to 2 r. Since for the earth, to is small 
(approximately 7.27 x 10' 5 rad/sec), the latter terms are most often neglected. The former 
is often retained for accuracy when computing the trajectory of a vehicle with near orbital 
capability. What is the maximum error that can be introduced into our analysis by 
neglecting either of these terms? For a given radial distance r, the transport acceleration 
depends on the latitude of the vehicle.The acceleration is zero at the poles and a maximum 
of value co2r when the vehicle is in the equatorial plane. The magnitude of this quantity is of 
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,he order 10-3 &> where g 0 is .he acceleration due to gravi.y a. .he Earth s surface. T e 
Coriolis acceleration depends on the magnitude and the direction of the vehicle's velocity 
with respect to the Earth. It is zero when the night path is parallel to the polar axis and a 
maximum when V is perpendicular to this axis. The magnitude of th.s value ts g.ven by 
2o)V and is of the order 10-> g„ at orbital velocities. Since we are primarily concerned wt. 
night in the sensible atmosphere for which velocities must remain sub-orbital due to 
aerodynamic heating constraints, we shall set CO equal to zero in the above equanons. Thts 

is the assumption of a non-rotating Earth [46]* 


Flat Earth. Further simplification results if the non-rotating Earth model discussed 
above is assumed fiat. Ignoring the term Wgr in the equations of motion amounts to 
ignoring the centrifugal force that contributes to the lift vector. As orbital veloc.ty ts 
approached the lift required is reduced to zero. Consequently the induced drag ts also 
reduced to zero. This is an important effect which must be accounted for. Thus, we shall 
not make the assumption of a flat Earth. 


Constant Acceleration due to Gravity. For flight within the atmosphere, the 
altitude of the vehicle remains small in comparison to the radius of the Earth. We may 
express the variation of acceleration due to gravity with altitude as g = g 0 r 0 2 / r where r - 
r + h. If we neglect h as small compared to r 0 , then g remains constant with increasing 
altitude. For an altitude of 300,000 feet, h / r 0 = 1.435 x 10‘ 2 and thus we induce an error 
less than that of neglecting the Coriolis acceleration. However, since including the variation 
of g with altitude is a trivial matter, the gravity field shall be taken as proportional to the 
inverse of the radial distance, r, squared. 


Constant Mass. This approximation is clearly invalid for a vehicle designed to 
achieve orbit. 1. is easily estimated assuming SCRAMJET propulsion and using an average 
specific impulse of 1500 seconds that half or more of the gross take-off we.ght for a 
vehicle with even a modest payload capability must be fuel weight (50J. 

Flight Constrained to a Vertical Plane. We shall find this approximation (bank 
angle of zero) to be most useful in getting started. However, it is expected that SCRAMJET 
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performance w„. require -ha- -he -rajec-ory fo.iow -he line of 
pressure, ahou- ». psf. Th-s is ,0 ensure a 
a result the vehicle will be operaong near minimum drag an 

;;:L, ve,oci, y is approached, iess and iess Hf. wi„ he required ” m"o 
Thrs may P-xxiuce -he need ,o fly a, negative hf, coefficM -o P=„ ^ ^ 
remain in the desired Mach number-alntude comdo M- 

consider abort and plane change maneuvers later in th,s analysrs. In e.-her case w 

from those of longitude, latitude, and heading. 

Angle of Attack and/or Flight Path Angle Small. I. is expected, for the 
reasons cited above in the discussion of flight constrained to a 
vehicleangleof attack - .remain 

srrjssim—— — 

will have to be evaluated numencally. 

Thrust Aligned with the Velocity Vector. This will be assumed until such time 
-hat a model for the normal component of SCRAMJET Unas, is available. A, very h,g 
Mach numbers, a large percentage of .he vehicle afterbody will ta«d io «pand 
exhaust now, reduce pressure drag, and increase the gross thrust, Stnce the pressu « ™ 
be" am above the free stieam pressure, they could add lift and increase the Ttft-drag 
ratio and consequently, cruise performance. In addition to making a large conmbut.on 
vehicle performance, propulsive lift may produce a nose-down p.tcta-g moment Tjn mus 
be trimmed ou, causing an aerodynamic penalty of increased mm drag. These fee 
prove .0 be very important and should be considered in later analys.s of th.s problem (39). 

Introduction of Specific Energy as a State Variable. This approximation 
w.ll be employed from the star, and offers order reduction through .be definition of a new 
Variable U 1.511- Its va.idity is iested in comparing the rescuing approx, mate soiu-.on 
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lo exact numerical results. Experience in actual nights and comparison between vanous 
oiut 1 or fighter aitcrafi have shown tha, the improvement in performance ,s — 
wtele exact optima, solemn is compttred with suboptima, soiuttons obtatned ustng thts 

approximation. 


Simplified Dynamic Model 

Implementing the previous, y cited assumptions numbered ,.5.8, and 9 in t^uanons (3- 

valid for aunospheric High, in a venical plane of a point mass over a sphertcal non-rot 
earth. The equations are: 


r = V siny 

(14) 

V(T - D) 
E = — : 

(15) 

m 


• i . u cosy V cos y 

7 " mV Vr 2 r 

(16) 

m = -f (r,V,T\) 

(17) 


It is assumed that the atmosphere is stationary and ^astate 1 ^^^ 

path (i.e. e t = 0). In equation (15), mass specific energy, E, is employed 


in place of velocity, V , where 
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E = 



l£ 

r 


hence V is to be taken as 


V = [2 (E + |i/r)] ,/2 


(18) 


(19) 


wherever it appears in this analysis unless otherwise noted. A dual-mode propulsion 
system is assumed. This system consists of a bank of SCRAMJET engine modules 
assumed to operate continuously at stoichiometric conditions above Mach 10 and a rocket 
engine that can be turned on or off as desired. The total fuel flow rate (f) and T are 
represented by 


T = T s + T) T r ; T] e [0,1] (20) 

f = c s T s + T) c r T r (21) 


where thrust specific fuel consumption is represented by c s for the SCRAMJET and c r for 
the rocket. The control variables are now rocket engine throttle (Tp and vehicle lift (L). The 
objective shall be to minimize the fuel consumed in an unconstrained energy climb. The 
performance index is given by, 

J = - m (tf ) ; tffree (22) 


5.2 Application of Singular Perturbation Theory 

The task is to determine the controls, T| and L, so that equation (22) is minimized. 
Even though the order of the dynamic model has been reduced from seven to four, a 
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TPBVP, that proves difficult to solve, again results from applying the necessary conditions 
for optimality. As stated earlier, this approach is computationally intensive and not well 
suited to on-board real-time implementation. For this reason we seek further model order 
reduction through the application of singular perturbation theory. This technique is detailed 
in references [12-14], Regarding energy and mass as slow variables and altitude and flight 
path angle as fast, we introduce the perturbation parameter e, nominally one, in equations 
(14-17). 


E = 


V(T — D) 


m 


m = -f (r,V,T|) 


’ L U cosy V cos y 
ey = — - ----- - + 


mV 


Vr 


er = V siny 


The first-order necessary conditions for optimality are: 


H = X e E + A. m m + A. r Vsiny + A^y = 0 


B~ 


an 

dE 


K=-~ 

dm 




(23) 

(24) 

(25) 

(26) 

(27) 

(28) 

(29) 




3H 

0r 


( 30 ) 
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(31) 


The analysis proceeds by examining the necessary conditions for optimality on two 
separate time scales, t and x = t/E, for £ = 0. The results are termed the reduced and 

boundary layer solutions, respectively. 

5.2.1 Reduced Solution 


The Unconstrained Case 


Setting e = 0 in (23-26) reduces the order of the dynamic system to two. Altitude and 
flight path angle dynamics are assumed fast in comparison to energy and mass dynamics, 
and altitude now takes on the role of a control variable. The differential equations (25) and 
(26) are reduced to algebraic equations which yield the following relations: 


70 =° 


L 0 =m 


£_v_ 

2 r 


(32) 


(33) 


The subscripted zeros denote reduced solution values in this context and are omitted where 
not deemed necessary for clarity. Thus, in the reduced solution, flight path angle is 
assumed to remain zero along the optimal path and lift is completely determined as weight 
minus centrifugal force. The reduced solution Hamiltonian is given by 



H o =A. E E + X m rn = 0 

(34) 

where 

O 

1 

II 

U-. 

' — * 

J 

(35) 


We can use the fact that H 0 is a constant of motion to express one of the unknown costates 
in terms of the other as follows. 



45 




f m 


v(t-d) 


(36) 


Note that when considering an energy climb (i.e. a positive energy rale) Xe must always 

have the same sign as V Furihetmore stnce ». represen. the variation™ 

index J with respect to mass, it cannot change sign (i.e. it is not possible for a reduction 

vehicle mass due to fuel expenditure along the climb pad, to incre^e 
vehicle). Given that J = -m(t f ). X„, and thus Xe. must rematn negattve. Using this i sign 

information it can be shown (see Appendix B) tha, maximizatton of the ' J 

divided by dm/dt with respect to the controls is equivalent to m.mmtzing the reduced 
solution Hamiltonian, (34), with respect to the controls. Thus, sattsfactton o e mmtmum 
principle is reduced to the following operation. 


V(T-D) 


h 0 > T lo =ar g max M f 


(37) 


E = constant 


In addition to being relieved of having to solve a TPBVP, the two-dimensional se^ch 
for values of rocket throtde setting and altitude that maximize the functton ' dennfle 
equation (37) while holding energy constant and constraining thrust to remain grea e 
drag can be performed off-line and stored as a function of energy level and mass. The u 
of equation (33) in such an off-line calculation requires that an estimate of vehicle a* T* 
generated as a function of energy level. A means fordoing sois reported on page 44 of dus 

report It turns out that with the current vehicle model being employed, e re “ c 
Itively insensitive ,0 mass variation. This is due to the fact that the vehtc e tends to 
operate a, very low lift coefftcients. As a result, when a change in mass effects the required 
lift through equation (33), the change in drag, which effects the maxim, zatton in (37 . 

only slight since the lift curve slope with respect to drag is very shallow at low values 

coefficient (see Figure A3). 


Note that rocket throttle setting, t). appears linearly in the Hamiltonian. This indicates a 
bang-bang control solution for rocket throttle setting and the possibility that stngular arcs 
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exist . I. is shown in the next sub-section of this report that intermediate values of roc et 
throttle setting, as would occur along an optimizing singular arc, are in fact non-optima . 
Furthermore, in the sub-section .ha, follows, an analytic rocket thro.de sw.tch.ng funcuon 
is derived that eliminates the need to perform the maximization tndtca.ed in equa 
with resnect to throttle setting. Thus, equation (37) reduces to the form: 


vcr-D) 

h 0 = arg max h ^ 


J E = constant 

T> D 

T1=T1* 


(38) 


Examination of Possible Singular Arcs 


Because thro.de setring appears linearly in rite Handhonian, the necessary contrition that 
the partial derivative of H, with respect to dtro.de setting be zero ytelds no control so u . 
directly. Instead we must examine the sign of this partial derivative m order to d * er ™ 
the optimal dtro.de setting. The optimal control necessarily occurs a, one control boundmy 
or the other and may switch from one to the other any number of times. There ts also . e 
possibility that singular arcs exists. That is, we may find that the afore mentioned P™ 
derivative is zero over a finite duration of time, hence no sign can be detemuned. In sue 
case intermediate values of throtde setting may occur and it must be determine w e 
such values are optimal. The switching condition (S) is determined as follows 


3Ho 

dn 



where T r > 0 


(39) 


The sign of the partial derivative of Ho with respect to Ti is determined by the sign of the 
last bracketed term in equation (39), which is termed the switching function, S. 


dH 


3ti 


S = X E (^)-X m c r 


(40) 
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The rocket control solution can be summarized as follows: 

T| = 1 when S < 0 

0>ti< 1 if S = 0 for a finite duration of time (41) 

t| = 0 when S > 0 

Next we would like to show that singular arcs, although they do exist in this problem, 
are not optimizing. By definition, the case in which an extremal arc (H u = 0) occurs along 
which the matrix H uu is singular is termed a singular arc. Such arcs satisfy the necessary 
condition on convexity, but not the strengthened condition; that is, H uu is only positive 
semi -definite. Additional tests are usually required to determine if a singular arc is 
optimizing or not [52,53], In this case it can simply be shown that H uu is not positive 
semidefinite. Thus, singular arcs are non-optimizing. We proceed as follows. 


Consider again the reduced solution Hamiltonian: 


H o =M + X m m = 0 


(42) 


Let us examine H uu where u T = [ V,q ]. Note that in equations (14-17), velocity, V, was 
eliminated through equation (19) when implementing the energy state approximation. An 
equivalent formulation of the reduced solution reported above is obtained by eliminating 
radial distance, r, and retaining velocity, V. In that case, V rather than h, acts as a control- 
like variable in the reduced solution. It is much more convenient to use the latter 
formulation for this argument. In this context, the subscript notation, H UUi denotes partial 
differentiation as follows: 


H 


*y 


. . T 

_a/9H\ 
3x\9y ] 


(43) 


Expanding H uu into its matrix representation yields. 


H 


111 


H w 
Hv n 


H„ v 


(44) 
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Recall that it is necessary for optimality that H uu be at least positive semidefinite. The 
determinant of the symmetric portion of H uu must be greater than or equal to zero for 
positive semidefiniteness. In the case of a singular arc, = 0 by definition. Thus, the 
only terms that contribute to the determinant of the symmetric portion of H uu are the off- 
diagonal terms of the matrix shown in equation (44). 


det H. 


Symmetric 


( H Tlvf + 


2 H^yHy-n 


(Hvt] 


(45) 


Carrying out the indicated partial differentiation, and ignoring the weak dependence of T r 
and C r on h yields: 


It is straightforward 


dH 0 


2 

aH 0 _ 

X E T r 

aqav 

m 

to show that 


2 

2 

a H 0 _ 

a h d 


av an an av 


1 c.Tj 

m r T 


(46) 


(47) 


Equation (47) also holds if the dependence of T r and C r on h is accounted for. Thus 

detH^ .- = - { 4 Hvn) 

111 symmetric ' 

which is negative for non-zero values of H Vtv The expression for H Vt1 is as follows 
where the approximation results from neglecting the weak dependence of T r and C r on h 
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Hv„= 


(v(t-d) Cfv ) Tr + ((t-d) 


_Cr T ry 


VTr 

m 


(49) 


In general the bracketed term of (49) is non zero. It follows from the condition H 0 = 0 that 
if Xe = 0 then X. m = 0, and visa versa, which is a trivial case. Thus H uu is not positive 
semidefinite, and intermediate values of rocket throttle setting are not optimal [54,55]. 


Analytic Form for Throttle Switching Condition 

Consider again the switching condition given in equation (40), 

^ = (“) - ^m c r (50) 


Using equation (36) we may eliminate the costate from equation (50) to yield. 


S =X 


m 


f 


T-D 



(51) 


As stated previously, X. m represents the variation in the performance index, J, with respect 
to mass, and X. m (tf) = - 1. Clearly, since J = - m(tf) , cannot change sign (i.e. it is not 
possible for a reduction in vehicle mass as fuel is expended along the climb path to increase 
the final mass of the vehicle). Therefore X. m must always be less than zero. The sign of Htj 
is thus determined by the bracketed term in equation (51). Substituting for thrust, T, and 
fuel flow rate, f, using equations (20) and (21), and taking into account (41), yields the 
following analytic form for the rocket throttle switching condition. 
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T| =0 



> D 


n = i 



< D 


(52) 


Estimation of States and Costates in the Reduced Solution 

In the above solution, the variables L 0 and y 0 are first eliminated before the 
maximization in (38) is performed. An alternative viewpoint is to adjoin the constraints 
(that the nght hand sides of (25) and (26) are zero) to the Hamiltonian. That is 

H 0 = X e 0 e + + \Vsiny + X y y (53) 


Associated with this formulation are the additional necessary conditions (see 30 and 31): 



II 

;ik 

8H °=0 

3L 

(54) 

which result in 

X =0 
X E 0 

2KV 2 L Q 

. q s 


(55) 

(56) 


Using (34) we may solve for X ^ in terms of Xm 0 



fm 

v(t-d) 


(57) 


where X mo satisfies (see equation 28) 
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^n=~ 


3Ho 

5m 


= X, 


V_\ 

2 

\m / 


T- D + - 


2KL ( 

qs 


(58) 


An approximate integral of (58) can be obtained by noting that the term KL D 2 /qs in (58) 
represents the induced drag component (see equation 12), which for small angles of attack 
is small in comparison to the zero lift drag component. Ignoring this term and using (57) to 
eliminate Xe, (58) becomes 



where again as in (35) Xm(tf) = — 1.0 

m 


(59) 


Integrating both sides leads to the approximation 



m(t f ) 

m 


(60) 


Note that expressions (56), (57), and (60) each depend on knowledge of m(tf). An 
estimate of m(t f ) can be formed by integration of dm/dt over dE/dt along the reduced 
solution as follows: 



( 61 ) 


then 


m(tf)= m(t 0 ) - Am 


(62) 


In a similar fashion it is possible to estimate the time required to gain energy. 


AT = 



dE 



dE 

E 


(63) 
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It is also of interest to form an estimate of the flight path angle along the reduced solution. 
This may be accomplished by first solving for the flight path angle, y, in equation (14). 


h = Vsiny 



(64) 


Using the chain rule for differentiation one can write. 


h = Jp-E 
dE 


(65) 


where dh/dE can be approximated by considering two closely spaced energy levels 


dh 

dE 


62 ' hi 
AE 


( 66 ) 


In equation (66) the superscripted asterisks denotes the reduced solution optimal altitude for 
a given energy level. Using equations (65) and (66) in (64), we find that the flight path 
angle at a given energy level, E 2 , can be estimated as follows. 


?2 = 


sin 


-1 


h 2 - hi E 2 
AE V 


(67) 


Inclusion of a Maximum Dynamic Pressure Constraint 

Recall the definition given for dynamic pressure, q, and velocity, V: 

q = l/2pV 2 = q(E,r) where V = [2(E + |i/r)] 1/2 (68) 

We wish to constrain q so that q < q m ax- Thus we wish to enforce a constraint on a 
function of the state and control variables (E and r [i.e., h], respectively) all along the 
reduced solution. That is 
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Minimize J, subject to the constraint, C (E,h) = q - q max < 0 (69) 

In order to proceed, we adjoin this constraint to the reduced solution Hamiltonian with a 
Lagrange multiplier as follows [52], 


where 


H o =X E E + X m m + A. q C = 0 ;X m (t f ) = -1.0 (70) 


X q =0, C<0 
x q >o, C = 0 


( 71 ) 


This constraint does not affect the throttle switching condition in (52), but it does limit the 
search space over which the maximization of equation (38) takes place: 


. V(T-D) 

h 0 ~ arg max h - 


£ = constant 
T> D 
T1=T1* 
q < qmax 


(72) 


The constraint multiplier Xq can be evaluated from the following condition. 


dH 0 

3h 


= 0 


9H 0 



ah 


(73) 


Reduced Solution: Numerical Results 

Reduced solution trajectories for the vehicle modeled in Sections 2-4 were generated by 
carrying out the maximization process indicated in equations (38) and (72) over the energy 
range corresponding to V = 5000 ft/sec at h = 0 to V = 25,000 ft/sec at h = 200,000 ft. The 
unconstrained case and the case for which dynamic pressure is constrained to be less than 
2000 psf are presented in Figure 21 . 



The unconstrained case begins at an altitude of 36,000 feet and initially exhibits an 
acceleration at constant altitude. This is followed by a rapid transition in altitude at 
approximately Mach 10 due in part to the modeling of atmospheric properties (an 
isothermal layer is assumed from h = 36,000 ft to 82,000 ft) and in part to a transition in 
SCRAMJET performance. Below Mach 10 optimum engine performance is achieved at a 
fuel-to-air ratio less than stoichiometric whereas above Mach 10 the calculated fuel-to-air 
ratio exceeds its stoichiometric value. It is assumed that no additional thrust is available at a 
higher ratio, thus the fuel-to-air ratio is constrained in the engine model to remain less than 
or equal to its stoichiometric value. This is indicated graphically in Figure 10 on page 22. 
From this point in the trajectory a slow steady climb ensues until the altitude exceeds 
82,000 ft, at which point a constant altitude acceleration is again optimal (since the assumed 
isothermal layer terminates). Note that at a velocity of approximately 21,000 ft/sec the 
rocket engine is deemed advantageous and is turned on. 

At a velocity of 22,000 ft/sec. operational efficiency dictates a switch in operating 
altitude. This is indicated in Figure 21 by a vertical line, and represents a constant energy 
instantaneous altitude transition. This altitude discontinuity also occurs in the q constrained 
case and in the case of SCRAMJET propulsion alone (no rocket) as illustrated in Figure 22. 
The mechanism causing this discontinuity is depicted in Figure 23 where the function to be 
maximized with respect to altitude[see equations (38) and (72)] is plotted as a function of 
altitude for three particular energy levels. Note that these curves terminate on the left at 
approximately 135,000 feet - this corresponds to the dynamic pressure constraint 
boundary. As the energy level increases, the local maximum at approximately 200,000 feet 
overcomes the functional maximum previously occurring at the constraint boundary. The 
shape of the curves presented in Figure 23 results from the interplay of the modeled thrust, 
drag, and fuel flow dependencies on altitude and velocity. The addition of rocket thrust 
shifts this altitude discontinuity to a lower energy level. The basic result is that at a 
sufficiently high Mach number, degrading engine performance and increasing drag causes a 
switch in the optimum operating conditions from high thrust/high drag to low thrust/low 
drag. Given an initial weight of 200,000 lbs, the weight of fuel expended along the 
unconstrained path shown in Figure 21 is estimated as 115,875 lbs. (58% of the take-off 
gross weight). The time of flight was estimated as 659 seconds. The most efficient cruise 
point was identified at approximately Mach 10 at 36,000 feet. The percent of take-off gross 
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r r ,vn» <;rR AMJET only case was estimated 

rrr. 

Optimal 

mass flow of air is available. Co . , 2000 psf Maximization 

^t“r,r:=r.rrsr:^*-<-”- 

of (72) subject to q < 2000 psf results m V - 21 500 ft/sec, (the rocket 

Here the vehicle follows the q consffaint all ^ which^point climb out of the 

being advantageous tiov' ^ cjted for „„ co „s.rained case. Note that 

atmosphere is initiated fo enforced! may force climb out 

temperature limits on the nose and wing leadmg edges > 

for the same — of the total gress weigh, being expended in achieving 
This corresponds to 2 percent mo the time 

th e specified final energy over the — ned - * alo „g the 

required to achieve the final energy, est.ma.ed as 37 • unco n Sff ained case is 

constrained path. The difference m time between tec accc ierate the 

;; h p — , C y Z“ but now on the q constrain, herniary at about 95,000 fee,. 

Pigures 24 and 25 presen. the required lift ^"lt 

attack required to mam, tun ''^“^'^"“e vehicle operates a, very low values of 
case as a function of energy • . , than f our degrees over the 

lift coefficient. The ^ ° .^.hras. to drag ratios computed along the 
entire trajectory. Ftgures 26 P , hat ^ average thrust to drag 

unconstrained and q constiamed cases Of nore^ ^ ^ ^ to ^ ratio for 

rano of the entire energy c ange fligh , path angle history required to 

the SCRAMJET only case is 1.62. The estimate g V ^ ^ spike a , 

follow the q constrained reduced solution chm pa ,s p earlier and even in this 

energy level 36 is genera, ed by the altitude discon, inu.t discus e ear ^ 

region, the High, path angle is estimated to for reliuced 

assumption of a small flight path angle and small angle of attack 

solution climb path. 
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Figure 21 . Reduced Solution, Unconstrained and q Constrained Fuel-Optimal Climb Paths 










Function being Maximized 


57 



Velocity in ft/sec 

Figure 22. Reduced Solution, q constrained fuel-optimal ascent 
trajectories. 
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Figure 23. Reduced Solution, mechanism by which altitude discontinuity 
occurs for the q constrained case, SCRAMJET only. 
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Energy Level 

Figure 24. Reduced Solution, required lift coefficient computed along 
the q constrained reduced solution. 
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Figure 25. Reduced Solution, required angle of attack computed along 
the q constrained reduced solution. 
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Energy Level 

Figure 26. Reduced Solution, Thrust to Drag ratio computed along the 
unconstrained reduced solution. ° * 




Energy Level 

Figure 27. Reduced Solution, Thrust to Drag ratio computed along the 
q constrained reduced solution. 
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Figure 28. Reduced Solution, History of Flight Path Angle estimated 
along q constrained path, SCRAMJET/rocket case & 


5.2.2 Boundary Layer Analysis 


The Unconstrained Case 

The unconstrained boundary layer solution associated with equations (23-26) is 
obtained by introducing the dme transformation x = t/e and again setting e = 0. The 

necessary conditions for optimality [56] become (where the prime notation denotes 
differentiation with respect to x ): 

h bl = ^e o E’ + X^m' + X r Vsiny + X y y' = 0 (74) 

E = 0 => E(t) = E 0 X’ E = 0 => \ e (t) = X Eo (75) 
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m' = 0 => m(x) = rrio X’ m = 0 => ^(t) = X n 


V- 


3H 


BL 


dy 


V s - 


3H 


BL 


3r 


r’ = V siny 


, L |i cosy V cosy 


y = 


mV Vr 2 


du 


BL 


= 0 


9L 


(76) 

(77) 

(78) 

(79) 

(80) 


Evaluation of the partial derivative in (80) results in the following expression. 


d H BL 

dL 



(81) 


Substituting equation (81) into equation (80) and solving for L yields the following 
expression for optimal lift (L*), 



The control solution for rocket throttle setting, r\, remains the same as in the reduced 
solution (see equations 52). 

The evaluation of Xy needed in equation (82) unfortunately requires the solution of a 
two-point boundary-value problem (TPBVP). The computational load associated with a 
TPBVP is to be avoided, and as a result we seek approximations to Xythat yield an 
acceptable guidance law. 



62 


SP1 Control Solution 

When close ,0 the reduced solution, it may be possible to use the expression for X* 
given by equation (56) as an approximation to \ Substitution of this zerot -or er 

approximation into (82) results in the following simple feedback law, 

L* = (p/p 0 ) U <83) 


where L c is given by equation (33) and p 0 represents atmospheric density at h h„. Note 
that this solution is independent of the reduced solution costa.es and the final mass. For 
convenience this control law is referred to as the SP1 solution in the figures that follow A 
numerical simulation using this control law to mack a constant energy condition reveals, 
however, that this solution is unstable. This is clearly illustrated in Figures and 
where an altitude of 36,000 feet and zero flight path angle are to be tracked given an mitia 
perturbation in altitude of 14,000 feet. This instability occurs for much smaller 
perturbations as well. Clearly we must seek a mom suitable approximation fork,. 


SP2 Control Solution 


Expansion of the boundary layer necessity conditions to first order about the reduced 
solution, which is an equilibrium point of these equations, results in die following linear 

perturbation equations [57-60]: 


’ 5r' 
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1 ** 1 
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( 84 ) 
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Time in seconds 

Figure 29. SP1 Controller Tracking Constant Energy Reduced Solution, 
Altitude versus Time 



Time in seconds 

Figure 30. SP1 Controller Tracking Constant Energy Reduced Solution, 
Flight Path Angle versus Time 
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The square matrix in equation (84) is a 
readily be evaluated analytically. 


A 2 j - 


Hamiltonian matrix, and all of its terms but A 31 can 


(85) 
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( 88 ) 


where H'bl is defined as in equation (74), but after elimination of the lift control us.ng 
(82) The rodte. throttle control is completely determined by the switchmg condmon gtven 
in equation (52, and Urns may be treated as a known constant. The ***** £ 
equation (87) can be evaluated numerically as follows. Constder a Taylors senes 

expansion of (74) to second order 


, . ' dHsti ... sWi 5h 

H B L(l>o + 8h ) = H BL ! + “ 5h + r- + . • • 


(89) 


3h i -V, 2 ; 2 


Taking into account the fact that H bl 1S zer0 
notation we have 2 

Hrl = H h 5h + H hh Sh /2 


dh‘ 

O 

along the optimal path and simplifying the 

(90) 


Define forward and backward perturbations with respect to h about h„, denoted + and 
respectively, as 
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H + = H h 5h + H hh 8h/2 (91) 

H" = -H h 5h + H hh 5h 2 /2 (92) 


Adding together H+ and H“ and then solving for the desired quantity we arrive at the 
following numerical estimate 


= H hh = (H%H")/8h 2 (93) 

3r 2 i 

l 

O 

A well-known property of Hamiltonian Matrices is that their eigenvalues are arranged 
symmetrically about the imaginary axis. Since we require that the boundary layer dynamics 
be stable forward in time, the state vector in (84) is expressed as. 


where 


x = k t a + k 2 b 


(94). 



8r, 8y, 8A, r , 81^ 


(95) 


and a and b are the real and imaginary parts of the eigenvectors associated with the stable 
eigenvalues. At any given time the state perturbations away from the reduced solution, 8r 
and Sy, are known. The scalars ki and k 2 are unknown, as are the perturbations 8X r and 
8Xy. Thus, equation (94) constitutes a system of four scalar equations in four unknowns. 
Use of equation (94) to pick the free values of 8X r and Sky allows the suppression of the 
unstable modes. That is, we allow only initial conditions that lie in the subspace spanned 
by the eigenvectors associated with the stable eigenvalues. To solve for kj and k 2 , 
partition equation (94) into upper, (U), and lower, (L), systems as follows. 







x u 

— * 

“ ^1 

a u 

+ k 2 

t>U 

X L _ 

a L_ 


v 


(96) 


Now since Sr and 8y are known, we can solve the upper system for ki and k 2 as follows. 
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x u - 


5r 

5y 


= kj ajj+ k2b(j 


(97) 


or 



(98) 


With ki and k2 known, it is a simple matter to calculate 5X r and 5X, y using the lower 
system, 



Now L* can be evaluated using equation (82) where 


( 100 ) 


That is. 



( 101 ) 


For convenience, this control law will be referred to as SP2 in the figures that follow. 


For the vehicle model being considered, the stable eigenvalues of (84) at the starting 
energy level are (- 0.2573 ± j 0.2599). These eigenvalues corresponds to a natural 
frequency ((On), damped natural frequency (a>d)> bme period (T), damping ratio (Q an< i 2% 
setding time (T s ) of: 

On = .3657 rad/sec T = 24. 1 8 sec 

co d = .2599 rad/sec T s (2%) = 15.54 sec 


; = .704 


( 102 ) 
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Figures 31 and 32 present altitude and night path angle lime histories generated using 
the control law (SP2) presented as equation (101). Excellent behavior, as predicted by 
(102), is evident in tracking the same constant energy conditions as in the previous Figures 
27 and 28. Difficulties arise, however, in that very large lift values (i.e. proportional to the 
altitude perturbation, which at 14,000 fee. is fairly large) are commanded. If a ltft limit is 
enforced, the performance of the SP2 controller is severely degraded. This condition is 
depicted in Figures 33 and 34 which again present altitude and flight path angle time 
histories for the same initial conditions, bu, with the magnitude of the SP2 commanded lift 
constrained so that L/W < 2.0. The corresponding lift control time histones for t e 
unconstrained and control constrained cases are presented as Figures 35 and 36, 
respectively. For the given altitude perturbation, the lift limited SP2 controller generates a 
20,000 foot undershoot. Clearly the lift limited SP2 control law is only suitable for small 
perturbations about the reduced solution. Configured as a launch vehicle, it is indeed 
reasonable to assume that the aerospace plane would not deviate far from the reduce 
solution climb path during ascent to orbit. However, two major stumbling blocks remain. 
Number one is the requirement that the vehicle traverse the large intenor altitude 
discontinuity present in the unconstrained (and the dynamic pressure communed) reduced 
solution climb path (see Figure 21 on page 56). This discontinuity introduces a very large 
altitude perturbation and would have to be traversed using an alternate controller. It has 
been shown that such problems can be addressed by constructing interior boundary layers 
[61,62], but again we are left with having to solve a TOBVP on line. Introduction of an 
aerodynamic heating constraint may eliminate the aldtude discontinuity from the reduced 
solution, but will still introduce a rapid change in slope in the climb path a. the juncture of 
the path with the constraint boundary. 

The second major stumbling block is the requirement for mission abort capability. 
Initiation of an abort maneuver may radically change the reduced solution and as a result the 
vehicle state may initially lie far from the reduced solution path. Again the approach o 
linearizing the boundary layer necessary conditions about the teduced solution may produce 

an unsuitable controller. 

In review, the SP1 and SP2 lift control solutions presented above were attempts to form 
a suitable approximation to the costate V hat is needed in equation (82) in order to 
alleviate having to solve a TPBVP. Note that control dependence on the final mass. 
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Figure 31. SP2 Controller Tracking Constant Energy Reduced Solution, 
Altitude versus Time 



Figure 32. SP2 Controller Tracking Constant Energy Reduced Solution, 
Flight Path Angle versus Time 
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Figure 33. SP2 Controller Tracking Constant Energy Reduced Solution, 
Altitude versus Time, LAV < 2.0 
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Figure 34. SP2 Controller Tracking Constant Energy Reduced Solution, 
Fligth Path Angle versus Time, LAV < 2.0 
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Figure 35. L/W for SP2 Controller Tracking Constant Energy Reduced Solution 
(see corresponding Figures 31 and 32) 
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Figure 36. L/W for SP2 Controller Tracking Constant Energy Reduced Solution, 
(see corresponding Figures 33 and 34) 
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although weak, is introduced through X yo (see 56, 57 and 60) in both SP1 and SP2. Recall 
that 5^ is determined in the SP2 control solution by linearization of the boundary layer 
necessary conditions about the reduced solution and that it has been shown that the 
resulting controller, when lift limited, is suitable only for small perturbations. As noted in a 
later subsection, state constraints in the boundary layer further complicate this linearization 
technique. It is suspected that in order to avoid the excessive lift commands that plague this 
linearized approach on must capture the full nonlinearites of the problem. The forced 
separation of altitude and flight path angle dynamics and a suitable penalty on flight path 
angle in the performance index as in reference [40] may yield a more accurate controller. It 
is proposed in Section 7 that this avenue be explored. Another problem cited in using the 
initial boundary layer lift control solution as the control law all along the trajectory is the 
lack of ability on the part of the controller to anticipate rapid changes in altitude; it can only 
react to such changes once they are encountered. In the next subsection of this report still 
another approach to constructing a lift control solution, admittedly suboptimal, is employed 
and shown to alleviate some of the afore mentioned problems. 


NLT Control Solution 

As an alternative approach to tracking the reduced solution, a nonlinear transformation 
(NLT) technique is employed as follows [63-68]. 

Again consider the altitude and flight path angle dynamics: 

i = V siny (103 ^ 

• _L it cosy V cosy (104 ) 

mV Vr 2 

Note that we have system equations in block triangular form, that is of the form 

Xi = f (X l , Xj) ( 105 ) 

x 2 = g(xi , x 2 , u) 


(106) 
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where xj = r, x 2 = y, and u = L. Next we take successive total time derivatives of r until 
explicit dependence on the control appears. 


r 


1 2 , ,2 2 
|asin y ^ Lcosy _ |icos y V cos y 

2 + m 2 r 

r r 


(107) 


Since the control, L, appears in the second time derivative we define U, the peusdo control, 
as 


U = x (108) 

It is desired that U be determined as follows 

U = K p (r 0 - r) + K d (r 0 -r) (109) 

where r 0 denotes the reduced solution radius at the current energy level and the time 
derivative of r Q denotes the climb rate required to stay on the reduced solution as energy is 
gained. This climb rate can be estimated by defining an appropriate increment in energy, 
evaluating the reduced solution at this higher energy level and then estimating the required 
climb rate using a forwards difference. 


The inverse transformation is defined by solving for L in (108) using (107) and (109), 


L = 



2 2 
V cos y 

m 

r 

cosy 


( 110 ) 


This lift control solution is referred to as the NLT control solution in the figures that 
follow. Note that as r and y approach their reduced solution values equation (110) 
approaches the reduced solution value of lift given by equation (33). A block diagram 
depicting the conceptual implementation of the nonlinear transformation technique to yield 
the controller defined by equation (1 10) is presented in Figure 37. This is mathematically 
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equivilent to the linear system depicted in Figure 38 which is used to design the controller. 
The corresponding closed loop transfer function is 


G(s) = 


K d s + K p 
s + K d s + Kp 


( 111 ) 


where the gains K p and K d for the second order system can be written in terms of the 
damping ratio, and natural frequency, ov as 

K p = (H2) 

K d =2Cco n (H3) 

The performance of this controller can be dictated by selecting the values of K p and IQ to 
yield the desired dynamic response. 

Figures 39 and 40 present the altitude and flight path angle time histories generated 
using the NLT controller to track the same constant energy conditions as before: h 0 = 
36,000 feet, y 0 = 0. The gains were selected so as to match the performance of the SP2 
controller at the same energy level. As expected, the trajectories and the corresponding lift 
control (Figure 43) are nearly identical to those of the SP2 controller that were given in 
Figures 31, 32 and 35. Figures 41 and 42 present the altitude and flight path angle histories 
for the case in which the magnitude of the NLT lift control is constrained so that LAV < 
2.0, for the same gains and for the same initial conditions as for the SP2 controller. The 
corresponding lift control is shown in Figure 44. As in the unconstrained case, the NLT 
controller very nearly duplicates the behavior of the SP2 controller. 

Figures 45, 46 and 47 present the case in which the NLT controller gains are modified 
to yield a 2% settling time of 120.0 seconds rather than the 15.5 seconds demanded by 
SP2. The resulting trajectories are well behaved. More importantly, LAV remains positive 
and less than 1.2. The NLT controller clearly provides a simple means for adjusting its 
performance to yield a control solution within reasonable bounds. It also provides 
guaranteed stability properties which the SP2 soludon lacks given large perturbations. 
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Figure 38. Block Diagram of Equilent Linear System 
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Figures 31 and 32) 




Time in seconds 

Figure 41. NLT Tracking Constant Energy Reduced Solution, 
Altitude versus Time, L/W < 2.0, Gains selected to match SP2 
(see corresponding Figures 33 and 34) 



Time in seconds 

Figure 42. NLT Tracking Constant Energy Reduced Solution, 
Flight Path Angle versus Time, L/W < 2.0, Gains selected 
to match SP2 (see corresponding Figures 33 and 34) 
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Time in seconds 

Figure 43. LAV for NLT Controller Tracking Constant Energy Reduced Solution, 
Gains selected to match SP2 (see corresponding Figures 39 and 40) 
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Figure 44. NLT Controller Tracking Constant Energy Reduced Solution, 

LAV < 2.0, Gains selected to match SP2 (see corresponding Figures 41 and 42) 
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Time in seconds 

Figure 45. NLT Controller Tracking Constant Energy Reduced Solution, 
Altitude versus Time, Gains picked to give a damping ratio of 0.7 and a 2% 
settling time of 120.0 seconds 



Time in seconds 

Figure 46. NLT Controller Tracking Constant Energy Reduced Solution, 
Flight Path Angle versus Time, Gains picked to give a damping ratio 
of 0.7 and a 2% settling time of 120.0 seconds 
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Figure 47. L/W for NLT Controller Tracking Constant Energy Reduced 
Solution, Gains picked to give a damping ratio of 0.7 and a 2 % settling time 
of 120.0 seconds (see corresponding Figures 45 and 46) 


Inclusion of a Maximum Dynamic Pressure Constraint 

Recall the definition given for dynamic pressure, q, and velocity, V: 

q = l/2pV 2 = q(E c 4) where V = [2 (Eq + |i/r)] 1/2 (114) 

It is desired to constrain q in the boundary layer so that q < qmax- Note that in the 
boundary layer, E is held constant at its reduced solution value and r resumes it status as a 
state variable. Thus, in contrast to the reduced solution in which the dynamic pressure 
constraint was a function of state and control, in the boundary layer we must enforce a pure 
state constraint. The associated tangency conditions which must be met at a juncture 
between unconstrained and state constrained arcs, along with possible discontinuitites in 
the costates, the Hamiltonian, and the control [52, 69, 70], appear to make a straight 
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forward application of the available theory difficult at best. A means for handling a pure 
state constraint in the boundary layer is currently under investigation. 

The lift control solution derived via the nonlinear transformation technique discussed 
when considering the unconstrained case applies equally well to the dynamic pressure 
constrained case. Figures 48 and 49 present the altitude and flight path angle time histories 
generated using the NLT controller to track constant conditions associated with the dynamic 
pressure constrained case for several values of K p and K<j , as given in Table 2. The 
corresponding lift control histories are presented in Figure 50. As indicated in Figure 48, 
the damping ratio can be chosen to avoid violation of the constraint boundary. In the next 
section the NLT control law is used to numerically simulate ascent to orbit along the q 
constrained reduced solution fuel optimal climb path. 


Table 2 NLT Gain Variations Depicted in Figures 48-50 


c 

tOn (rad/sec) 

T s (sec)* 

K p 

Kd 

0.707 

0.047 

120.0 

0.0022 

0.0665 

0.854 

0.039 

120.0 

0.0015 

0.0666 

1.0 

0.033 

120.0 

0.0011 

0.0666 


* 2% Settling Time given by 4/^(0,, 



Flight Path Angle in degrees 


85000 





st 

*«•* 


<%> 

"a 

a 




80000 “ 

75000 “ 

70000 " 

65000 “ 

60000 L 
-50 



Time in seconds 


Figure 48. NLT Controller Tracking Constant Energy q Constrained 
Reduced Solution, Altitude versus Time, Gains picked to yield various 
damping ratios (see Table 2) but to maintain a 2% settling time of 120.0 
seconds 
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Figure 49. NLT Controller Tracking q Constrained Constant Energy 
Reduced Solution, Flight Path Angle versus Time, Gains selected to yield 
various damping ratios (see Table 2) but to maintain a 2% settling time of 
120.0 seconds 
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Figure 50. LAV for NLT Controller Tracking Constant Energy Reduced 
Solution, Gains selected to yield various damping ratios (seeTahle 2) but to 
maintain a 2% settling time of 120.0 seconds (see corresponding ig 
and 49) 
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6. Numerical Simulation Results 

A numerical simulation of the dynamic pressure constrained <q < 2000 
energy climb was earned out using the feedback law (NLT7 gtven tn equation ^ 1 0 wtd, 
gains selected to yield a damping ratio of .7 and a 2% setthng time of 1200 « 
lepurpose of comparison, thevalue of theft. 

iTgyleve. and the initial nigh, path angle was specified as zero. The resulting t^ctoty 
superimposed over the q constrained reduced solution, is presented along wtth the 
associated lift control time history in Figure 51. As expected, a slight lag ts ev, 
throughout the trajectory and a large value of lift is commanded upon reaching , e 
discontinuity. Figure 52 presents a second numerical simulation given the same tn.ttal 

conditions and controller gains but with the firs, time derivative of S 

forwards difference in energy equal to l/50th the final energy level. Excellent macta g 
th = reduced solution path is now achieved with the vehicle flying d**y 
dynamic pressure constrain, boundary. Note that in et.her case the add,, ton of rock 
propulsion is deemed advantageous via (52) and switched on. The flight path ang 

ny either of the trajectories shown in Ftgures 5. and 52 are presented as 

Figures 54 and 55, respectively. In either case i, remains less than ' 

assumption of zero flight path angle in the reduced solutton turns ou 
Approximately sixty percent of the total vehicle ma* is consumed tn ehmbtngto the 
energy given an initial weigh, of 200,000 pounds. This percentage of ^gross we tght 
consumed is equivalent to that estimated as required to follow the reduced solution ^ 
discussion on page 55) and i, thus appeam that funher optimization of the 

path angle dynamics is of no value in minimizing the performance index. g 

acceleration estimated as 0.6 g's when averaged over the entire « “ 

g’s through the Mach range M = 9 ,0 13 where the excess thrus, avertable ,s greatest. 

There is still a lag in the response a, the energy level where the reduced solution altitude 
is discontinuous. This lag can be avoided by calculating r„ and the desrred rate of cl.mt 
using a -forward look" procedure. This simply consists of performing tire calculauons of r„ 
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and dr^/dt at an energy level slightly higher (ahead of) the current energy leveh It should be 
possible to derive some intelligent logic .ha, adjusts the controller gatns and the magu tude 
-forward look" as a function of the current vehicle state. An tnvest.gatton of thts 

matter is proposed as a future research objective in the next section. 
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Time in seconds 

Figure 53. NLT Guided Solution, Flight Path Angle Time History, 
No estimate of desired rate of climb. 



Time in seconds 

Figure 54. NLT Guided Solution, Flight Path Angle Time History, 
with estimate of desired rate of climb. 
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7. Conclusions and Topics for Future Research 

This report documents the work accomplished during the period July 1987 to 
September 1988. All of the work documented in the progress report previously submitted 
for the period July 1987 to December 1987 is included in this report with several very 
important corrections implemented within the SCRAMJET performance code. These 
corrections removed several disparities reported between our results and those of cited 
references as follows: 

• The SCRAMJET model, based in large part on empirical relations, now exhibits 
qualitative behavior similar to that reported by other sources. 

• Sizing of the engine inlet area by estimating the vehicle's level flight envelope 
now indicates an increased requirement for inlet area as Mach number increases. 
This is the expected trend that was not observed previously due to incorrect 
SCRAMJET modeling. 

It was reported previously that the analytic aerodynamic vehicle model being developed 
was not yet a satisfactory tool for use in trajectory optimization and guidance law 
development. In particular it was stated that the hypersonic aerodynamic model 
satisfactorily predicts drag variation with lift and is responsive to control surface deflections 
but poorly estimates vehicle lift and its variation with angle of attack. The work completed 
in this area as of December 1987 is documented in Appendix A of this report . As cited in 
the appendix, suitable corrections to the method may take the form of additional geometric 
complexity but such have not been implemented at this time. 

A number of important modeling issues have not yet been addressed and many of these 
may well have a sizable influence on vehicle design, vehicle performance, and the nature of 
optimal trajectories. These include sensitivity of SCRAMJET performance to angle of 
attack variation, the incorporation of means for accounting for fuel pre-heat due to its 
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circulation as a coolant prior to combustion, and the addition of fuel into the combustor in 
excess of the stoichiometric ratio for the purpose of engine structural cooling. Another 
important effect not yet accounted for is the possibility of a large component of 
SCRAMJET thrust normal to the body longitudinal axis. This force, which contributes to 
overall vehicle lift and which can cause a large nose down pitching moment, may have a 
significant impact on vehicle performance and should be properly modeled. The vehicle 
model's aerodynamic envelope must also be enlarged to incorporate subsonic and 
supersonic flight and alternate modes of propulsion for these lower speeds. Additional 
controls such as variable geometry inlets may also need to be modeled in the future. 

It is not yet clear whether or not intermediate values of SCRAMJET thrust are desirable 
during ascent or an abort maneuver. Nor is it clear how a throttle back would be achieved. 
Most likely one or more engines modules would be shut down and their intakes closed to 
reduce drag. This leads to an optimal control problem in which the control variable is 
constrained to take on only a finite number of discrete values. The minimum principle still 
holds in such a case but the necessary conditions used to formulate the control solution 
presented in this report no longer apply. New necessary conditions would have to be 
derived to proceed with the analysis. 

Several trajectory constraints have been ignored in the analysis thus far. These include 
acceleration limits, angle of attack limits, and most importantly, aerodynamic heating 
constraints. Of note is the fact that the incorporation of these constraints will likely remove 
the large altitude discontinuity observed in the reduced solution in the unconstrained and 
dynamic pressure constrained cases. 

Consideration of either dynamic pressure or aerodynamic heating constraints leads to a 
state constrained boundary layer problem. Great difficulty has been encountered in 
attempting to construct a linearized boundary layer solution when a state constraint is to be 
enforced since no theory has been developed for such a case. The theoretical aspects of this 
issue need to be addressed and, if possible, a suitable methodology for handling such 
problems should be developed. In addition, alternate approaches, such as the forced 
separation of altitude and flight path angle dynamics, need to be investigated and compared 
with the linearized approach. The suboptimal nonlinear transformation approach outlined in 
Section 5.2.2 shows great promise with regard to overcoming some of the shortcomings of 
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the linearized approach. An extension of this technique to include transformation of the 
performance index so that a linear quadratic cost function appears in the linear space would 
result in a very powerful optimal approach and should also be investigated. 


Proposed future research tasks are detailed below. 

• Further develop the vehicle model as follows: 

- Extend the aerodynamic model to the subsonic and supersonic regimes and 
define a suitable means for transitioning between each flight regime's data 

set. 

- Develop suitable turbojet and ramjet engine models. 

- Investigate additional SCRAMJET modeling issues including variation of 
engine performance with angle of attack, additional thrust derived from fuel 
preheat, additional fuel flow required for engine structural cooling at the 
higher Mach numbers, a possible requirement for variable geometry, and 
dual-mode (ramjet-SCRAMJET) operation. 

- Incorporate the Generic Hypersonic Aerodynamic Model Example 
(GHAME) as an alternate aerodynamic model. 

- Develop suitable analytic (to the extent possible) aerodynamic models that 
will allow trajectory studies to be conducted on configurations for which 
windtunnel data is not available. 


• Extend the analysis to include: 

- Flight in the Mach range 0-5 (This task will require consideration of optimal 
propulsion system transitions in the subsonic and supersonic flight 

regimes.) 
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- The component of thrust normal to the flight path 

. Variable SCRAMJET thrust (This .ask may requtre ^ deriva.mn of 
necessary conditions associated with discrete variauon of SCRAMJE 
thrust since it is likely that SCRAMJET thrust variation will be achieved y 
turning on or off individual engine modules.) 

. Aerodynamic heating, acceleration, angle of attack, and lift constraints 


- Orbital insertion end conditions 

- Three dimensional dynamics as may be requirad for lift modulation, orbital 
plane change, and abort maneuvers 

- A rotating oblate Earth 

- A nonstationary atmosphere 


. Condnue the derivation of a suitable guidance algorithm. 

- Address the theoretical issues associated with the inclusion of a state 
variable constraint in the boundary layer analysis. If possible devise a 
method for synthesizing a linearized boundary layer lift control soludon for 

the constrained case. 

- Consider the faced separation of altitude and flight path angle dynamics as 
an alternate apptoach to handling the boundary layer problem. 

. Consider the derivation of intelligent rules for selecting gains and esumadng 
the desired rate of climb in the NLT lift control solution. 
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- Consider the possibility of an exact nonlinear transformation of the 
boundary layer necessary conditions to a linear optimal control problem 
with a quadratic index of performance. 

- Extend the Singular Perturbation formulation to include the out-of-vertical- 
plane dynamics as an additional layer. 


• Evaluate the resulting real-time guidance algorithms in nonreal-time simulation 
studies. 

- Compare the guided solutions with exact numerical solutions obtained using 
a multiple shooting algorithm. 

- Examine the robustness of the guided solutions to variations in atmospheric 
conditions, off-nominal engine performance, and other modeling 
uncertainties. 


• Conduct sensitivity studies to examine the impact of vehicle sizing parameters 
on the nature of the optimal trajectories. 
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Appendix A 


Prediction of Lift and Drag in Hypersonic Flow 


The flow features that dominate the aerodynamic behavior of atmospheric flight vehicles 
vary greatly with the flow velocity and the thermodynamic properties of the fluid medium. 
This variation is perhaps best quantified by means of the Mach number. The methods that 
may be employed in estimating aerodynamic forces in turn depend in large part on which 
flow features are dominant. For this reason it is convienient to divide the total flight regime 
into four regions, namely the subsonic, transonic, supersonic and hypersonic flow 
regimes. Table Al, on the following page, provides a definition of each regime and rules of 
thumb as to their relative boundaries [71]. 


There are, of course, other important ways of classifying flowfields. For example, 
flows in which the effects of viscosity, thermal conduction and mass diffusion are 
important are called viscous flows. Of note is the fact that surface pressure distributions, as 
well as aerodynamic lift and moments on some bodies can be accurately obtained by means 
of the assumption of inviscid flow. Another common assumption is that the fluid medium 
is a continium. This assumption is violated only for very low density flows, which occur at 
very high altitudes (above 200,000 feet) [72], This appendix is dedicated to aerodynamic 
force predictions in the hypersonic flight regime. Methods suitable for developing subsonic 
and supersonic models of similar complexity are currently being investigated. 

There is no clear dividing line between supersonic and hypersonic flow and the often 
quoted boundary of Mach 5 is, in reality, only a rule of thumb. Hypersonic flow is 
formally defined as that regime where one or more of the following phenomena dominate 
the flow field [73]: 

♦ Thin Shock Layers 

♦ Entropy Layer 

♦ Viscous Interaction 

♦ High Temperature Flows 

♦ Low Density Flows 
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Table A 1 

Definition of Flow Regions 


Flow Region 

Definition 

Rule of Thumb 

Subsonic-incompressible 

Density is constant 

Moo < 0.3 

Subsonic-compressible 

M < 1 everywhere 

Moo < 0.8 

Transonic Mixed regions where M < 1 and M > 1 

0.8 < Moo < 1.2 

Supersonic 

M > 1 everywhere 

Moo > 1.2 

Hypersonic 

See text that follows 

Moo > 5.0 


As one would suspect, these effects - thin shock layers and hot chemically reacting 
gases — add great complexity to the analysis of supersonic flows. In fact, the solutions to 
such problems push the limits of current technology [73]. Modem hypersonic research is 
now dominated by the methods of computational fluid dynamics (CFD) and without such 
tools and the supercomputers on which they are exercised it is unlikely that the aerospace 
plane would ever be bom [38]. 

Fortunately for our purposes (trajectory optimization) we find that viscous interaction, 
high temperature and low density effects may be disregarded. The approximate methods 
developed for inviscid hypersonic flow will allow us to adequately estimate lift and wave 
drag coefficients [73]. A major simplification results from the "Mach Number 
Independence Principle " which is illustrated in Fig. Al. This figure, reproduced from 
reference [73] and generated using the oblique shock relations, indicates that although the 
pressure ratio, Pj/P2, is continuously decreasing, the pressure coefficient, Cp, is 

approximately constant above Mach 6 or 7. Thus the force and moment coefficients 
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obtained by integrating the pressure distribution over the body are also independent of 
Mach number in the hypersonic regime. 

Additional simplifications accrue from the application of impact methods to represent the 
pressure distributions over the body and aerodynamic surfaces. Newtonian flow theory, 
which is further detailed at the end of this appendix, does not assume a continuim, but 
rather models the flow as a stream of discrete particles. It is then postulated that the normal 
component of momentum of each particle is destroyed upon impact with a body immersed 
in the flow whereas the tangential component is assumed to remain unchanged. The force 
exerted on a flat plate by the presence of the flow is then easily computed by using the 
conservation laws. 

Since the wings, fins and control surfaces of our vehicle are thin, they can be 
approximated as flat plates of zero thickness. Simple Newtonian flow theory then provides 
a simple analytic means for predicting the pressure distribution on these surfaces as a 
fountain of angle of attack and independent of Mach numbers greater than about 5. Of 
note, however, is the fact that at hypersonic speeds, the flow about a flat plate with a sharp 
leading edge "sees" a blunt nosed body due to the very rapid buildup of a thick boundary 
layer (i.e. viscous interaction). This is often referred to as the leading edge problem. 
Perhaps more to the point is the fact that all practical hypersonic vehicles have blunt noses 
and leading edges to reduce aerodynamic heating. Newtonian flow theory does not account 
for the additional drag due to a blunted leading edge when approximating a wing as a flat 
plate or the blunted nose of a fuselage. Blast wave theory, which is also further detailed at 
the end of this appendix, provides a fairly simple means for correcting the former method 
for these important leading edge and nose effects [74]. 

In some cases a combination of Newtonian flow and blast wave theories yields a highly 
accurate model of the pressure distribution on a body at hypersonic speeds. Figure A2, 
reproduced from reference [73], compares the pressure coefficients obtained using 
combined blast wave/Newtonian theory with flight data for the space shuttle. Remarkable 
agreement is achieved. 

This technique has been applied to develop a simple three-dimensional aerodynamic 
representation of a slender hypersonic vehicle using a minimum number of geometric 
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parameters. The resulting model, valid only in the hypersonic flight regime, can be adjusted 
in a straightforward fashion to simulate the geometry of a wide variety of aerospace plane 
configurations. This method has met with limited success when applied to the selected 
hypersonic research vehicle configuration. A three view drawing of the vehicle 
configuration was presented in Fig. 2 (of the main body of this report). The side, 
planform, and fin profiles were fit with straight line segments. Ten such segments yield an 
acceptable representation. The effects of wing incidence, fin cant, and propulsion modules 



Figure A1 Mach Number Independence 
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Figure A2 Comparison of pressure coefficients obtained with combined 
blast wave/Newtonian theory (solid line) with flight data for the space 
shuttle (circles). Windward centerline. Mach number of 21.6, angle ot 
attack of 40°. 


were incorporated. Blast wave theory was applied to the body to estimate the nose effects 
on the pressure distribution. Newtonian impact theory was used to correct this pressure 
distribution for yaw and pitch orientation angles (angle of attack and sideslip angle), and 
was also used to represent the wave lift and drag on the aerodynamic surfaces. Provision 
was made for differential deflection of the elevons. No attempt was made to correct for 
interference effects. Viscous drag was estimated by application of a simple skin friction 
coefficient. Figures A3 and A4 show the predicted drag polar and variation of L/D ratio 
with Cl for the model described. Comparison with experimental results at Mach 6 [18], 
indicates reasonable agreement between theoretical and measured results. 
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Figure A5 shows an important lateral stability effect and its dependence on the fin cant 
angle. The two curves representing 0 and 10 degree cant angles were generated by the 
method described above. Notice that the yaw moment coefficient versus sideslip angle 
exhibits a zero slope at zero sideslip. This results in poor lateral stability characteristics. 
This difficulty is alleviated by applying a cant angle to the fins as illustrated. 

Figure A6 shows the predicted effects of eleven deflection on the longitudinal 
aerodynamic behavior as represented by the pitching moment coefficient. Cm. Clearly the 
forces generated by elevon deflections are well represented, but the overall trend in Cm 
does not compare well with measured results. This is due to the very poor performance of 
the model in predicting vehicle lift. Figure A7 presents the trend of predicted and measured 
lift coefficient with increasing angle of attack for neutral elevon settings. Examination of 
this figure reveals the nature of the problem. The Newtonian result for the lift coefficient of 
a flat plate at angle of attack (see the end of this appendix for a derivation of this result) is: 

Cl = 2 sin^a I cos a I 

This trigometric behavior is clearly indicated in the predicted results by the inflection at 
zero angle of attack. On the other hand the measured results reflect a near linear relation 
between lift coefficient and angle of attack. Thus, although the model reasonably predicts 
drag as a function of lift, as required in the study of vehicle dynamics where lift is a control 
variable, the poor prediction of lift as a function of angle of attack precludes using the 
model to determine vehicle trim conditions. In order to allow work in trajectory 
optimization to continue, curve fits to the experimental data have been constructed and are 
to be used until a suitable model can be constructed. It is anticipated that the introduction of 
further geometric detail will result in a satisfactory model. The effect of body thickness, 
which has thus far been ignored in lift calculation, will provide additional lift generation at 
very low angles of attack. This mechanism is clearly evident in Fig. A5 by noting that the 
behavior of the vehicle with zero fin cant is similar to a body of parallel sides or zero 
thickness whereas the behavior of the vehicle with fin toe-in is similar to a tapered body. 
The result of the taper is more nearly linear behavior with changing incidence angle. 


Combined Newtonian Flow and Blast Wave Theory 


Simple Newtonian [74] 

This method provides simple but quite accurate estimates of surface pressure, exact as 
Moo, the free stream Mach number, tends to infinity and k,the ratio of specific heats, tends 
to 1.0. 

Consider a body (in this case a flat plate) immersed in a fluid medium with free stream 
Mach number much greater than 1. We consider the fluid not as a continium, but as 
composed of discrete particles that do not interact with one another. This was Newton's 
flow model. Newton postulated that the normal component of momentum of each fluid 
particle was destroyed upon impact with a body immersed in the flow. The tangential 
component of velocity was assumed to remain unchanged. Figure A.8 illustrates this 
model. The force exerted on the flat plate by the presence of the flow can be computed by 
using the conservation laws. Newton used this flow model in an attempt to explain the drag 
of projectiles. At low speeds this model is very poor (it predicts that Cl is proportional to 
a 2 when really Cl is proportional to a), but at very high speeds (i.e. in the hypersonic flow 
regime) it is quite accurate. 

Using Newton's flow model, it can be shown that the pressure coefficient on a flat 
surface with M co » 1 is approximated by: 

Cp = 2sin2a (A.l) 

This relation is exact as Mootends to infinity and k tends to 1.0 which in reality is never the 

case. This approximation is good for high Mach numbers however, continually improving 
with Mach number till M > 10 or so. Note that for air k = 1.4, but for M** » 1, k tends to 

decrease towards 1.0 as desired due to increasing temperature and its affects such as 
ionization. It is assumed that the pressure exerted on the surface of the body is zero 
everywhere in the aerodynamic shadow (see Fig. A.8). 
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Ut us apply simple Newtonian theory to estimate the lift and drag of a flat plate at angle 

of attack tn hypersonic flight. Referring to Fig. A.9. the force per unit span on the plate is 
given by: 


F ( c ) p lower - (c) p upper (A.2) 

t.e. the net force on the plate (the prime denotes per unit span) is equal to the pressure 
deference between the upper and lower surfaces dmes multiplied by the surface ruea. 

Now, adding and subtracting P„ (pressure in the free stream) we can force this relation 
into a form in terms of the pressure coefficient, defined as: 

C p = ( P - Poo) / ( 1/2 Poo Voo 2 ). ( A . 3) 

We then have: 


F - [ (Pj - p^) . (p u . pj j c (A.4) 

where P u - P„ is approximately zero since the upper surface is in the aerodynamic 
shadow. 

Then 

C F -= F / ( 1/2 p„ c ) = (P, ■ P M ) / (1/2 p M V M 2) = Cp, = 2 sin2 a (A.5) 

where C F ' denotes a non-dimensional force coefficient. Resolving this force coefficient 
into lift and drag components we have: 


Cl' = 2 sin2 a Icos al 


(A. 6) 


Cd’ = 2 sin 2 a Isin al 


(A.7) 




Figure A8 Newton’s Row Model 



Figure A9 Components 
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Blast Wave Theory [74] 

We may consider bias, wave .heory as a variation of .he hypersonic equivalence 

I hich slates- "A steady 3-D hypersonic flow over a slender body is equ.va en 
pr.nc.ple which sta.es. A s y yp body we may approximate 

an unsteady flow in one less d.mens.on. For y 

.he now situation by assuming all of .he drag is concen.ra.ed a. .he > 

We may then reduce a 3-D blunt-nosed body problem to that of a 2-D areolar shock 

grows with time (see Fig. A. 11). 

Much work was done on the theory of circular and spherical blast waves in the 1940's ; - 
,950 s These results can be used .0 approximate the shock shape for blunt-nosed bod e 

in hypersonic How fields by equating %*££££ £ ££££. dete^inel 
this technique requires knowledge of Cj 

the oblique shock relations are used to estimate the surface pressure distnbuuons. 

The blast wave solutions are based upon the assumption that flow similarity exists, i.e. 
it is assumed that 

P(r) / P(R) = function (r / R) 

where r is the distance to a point of interest in the field and R is the distance Tmm tiu ***** 

to the shock Very close to the origin of an expansion we expect extrem P 

rrrrrr 

far back from the nose. 

We will now develop the blast wave solution for a blunt nosed c >’ hn ' icr (, _' “ 
axisymmetric body) After noting that this solution is only available for zero . angle o ari ack, 

a correction for angle of attack based on simple Newtonian will be 
hypersonic speeds, the fiow about a fiat plate with a sharp lead.ng edge es a b 
h hodv due to the very rapid buildup of a thick boundary layer (t.e. 
i„°, "action) This is often referred to as the leadingedge problem. We may find it necessary 
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employ bias, wave theory to estimate the additional drag due to a 
when approximating a wing as a flat plate. 


blunted leading edge 


Blunt-Nosed Cylinder at Zero Angle of Attac 

For an expanding shock from a finite release of energy in a gas at rest blast wave theory 

givesthendiusofdtes^as:^ ^ ^ (A.8) 


where f (k,n) = constant from numerical integration which 
k = ratio of specific heats 
n = characteristic number 
n = 0; planar case 
n = 1; cylindrical case 
n = 2; spherical case 


Now we make a transformation from time to space: 

t becomes x/V„ and R(l) becomes R<*> 

R(x) = f(k,n) [E/p ■.)<>»«» “ M:M> (A ' 9) 

This gives the shoe, shape as a function of the iongitudinai body coordinate x (see Fig. 
A. 3). 

Next we equate energy 


in the blast problem to drag in the aerodynamic problem. 


E = D = qooC(j[rcd 2 /4] (n-1) 


(A. 10) 


Here d is the diameter of the cylinder and q^ is dynamic pressure. Then 

R/d = f (k, 1 ) { [1/2 PooV oo2 C d fit d 2 /4)] / (poo d 4 ] } 1/4 CA.11) 

R/d = f c (k) c d fr 4 [m m where W is e iven as 0795 (AJ2) 
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Thus if we know the vahte of C d , we can ge, an estimate of the shock shape. We «n 
then use the oblique shock relations to ge. C p jus. behind the shock. Frnally we assume 
Cp shock ' s approximately equal to Cp body. 

Pe/P„ = l2k/(k+l)]M„2 sin 2p- [(k-lV(k + l)] (A.13) 


We can neglect the second term of (13) for strong 
lies back very near the body we may assume: 


shocks and since for Moo» 1 the shock 


sin P = tan p = (3 = dR/dx = d(R/d)/d(x/d) 
sin2p = { d/d(x/d) [Eqn. (A. 12)] }2 


So that 


(A. 14) 


(A. 15) 


Pc/Poo = g(k) Moo 2 [ C d ^/(x/d)] = Pbody/ p « (A * 16) 


where g(k) is a constant given as 0.067. Now by definition: 


C p = 2/(k Moo 2 ) [ P/Poo - 1] 


(A. 17) 


If we neglect the 1 as small compared to the pressure ratio then we may write: 

C p = 2 g(k) / k [ Cd 1 / 2 / (x/d)] IUL] (A. 1 8) 

where L is total body length. Taking k to be 1 .4 the quantity 2 g(k> / k becomes 0.096 and 


where L is total bo 
we have the result: 


C p = 0.096 C d 1/2 (L/x) (L/d) - 1 


(A. 19) 


.,tr nrlc nf yero We may correct for incidence 
Note that this results assumes an angle of attack ot zero. y 

angle using the simple Newtonian result derived for a flat plate: 

Cp = 0.096 C d 1/2 (Ux) ( L / d ) - 1 + 2 sin2 a 


(A. 20) 



Appendix B 


Minimization of a Hamiltonian Function with One Unknown 
Costate and One Costate of Known Sign 


This appendix generalizes Appendix E of reference [40] and as such documents a 
method for minimizing a Hamiltonian with one unknown costate and one costate of known 

sign. 


Consider the Hamiltonian function 

H*(x) = \i f(x) + X g(x) = 0 


(B.l) 


where the sign of p is known. It is of interest to minimize this function . with respect to *. 
Alternatively we may consider an equivalent operation exerctsed on the funcnon H(x) w, 
respect to x where H(x) is defined as. 


H(x) 


= H*(x)/ji = f(x) + n g(x) = 0 (B.2) 


and for convenience we have defined. 


T\ = X/\l 


(B.3) 


Now, equivalent to minimizing 


H*(x) with respect to x we have the following operations 


min* H 


if p > 0 


max* H 


ifp < 0 


(note p.^0) (B.4) 




B2 

— 

Sufficient conditions for the existence of a minimum are (where the subscript notation 


denotes partial differentiation) 



H x = f x +t| g x = 0 

(B.5) 

— 

f > 0 if H > 0 



Hxx = fxx + B gxx < 

(B.6) 


^<0 if p. < 0 


— 

In a free final time problem where H* does not explicitly depend on time, it is also 
necessary that H*. and thus H, be zero. Using (B.2) to solve for r\ leads to 

— 

n = -f/g g*o 

(B.7) 

— 

Using (B.7) in (B.5) and (B.6) we obtain the conditions: 


— 

H x = fx + (~f/g)g* = 0 

(B.8) 


V 
o 

»-*• 

*-n 

T= 

V 

0 


Hxx = fxx “*" ( f /§) gxx \ 

(B.9) 


[<0 if n < 

0 


Now define the function 


— 

L = g/f 

(B.10) 

— 

Next take the first and second partial derivatives of L with respect to x 



Lx = (f gx “ g f x) / f 2 

(B.ll) 


L X x = (f gxx “ g fxx) / f 2 ~ 2 f x L x / f 

(B.12) 

— 

Setting (B.ll) equal to zero yields the same conditions as given by (B.5) and (B.7). 
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Condition (B 9 ) ran h» 

K ) can be rewritten as 

gf « " f gxx >0 if g>0 ) 
gf «~ f gxx<0 if n 


^ingL x = 0in( B .i 2 ) 


g < 0 


we have 


>0 Sf "' f 8«<0 if g>0 

if g<0 j 


4 < 0 


(B.13) 


LXX ~ (gf « ~ f 8xx)(~l /f2) 


Since the last term in rR • , 

CB. 13): ' } 1S 3 Ways "native, the followin 


(B.14) 

g conditions are equivalent to 
Lxx > 0 if g > o) 


Lxx>° if g < oj 

i. <- o * t 4 > 0 

From the fore8 ™ & we ca " «*** to 
max x {L} if g > o ) 

nun * (LJ if g > 0 j) 


CB.15) 


“in* (L) if g< o M>0 


max * {L} if g<0 j 


f4<0 


(B.I6) 


C e< 3 u * Va Jcnt to the conditions CB 2 n * D 


S'g” or that (i and X always change 


sign 



